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ABSTRACT 



We report on Herschel/PACS observations of absorption lines of OH + , ILO+ and H 3 + in NGC 4418 and Arp 220. Excited lines 
of OH + and H 2 + with £i ower of at least 285 and ~ 200 K, respectively, are detected in both sources, indicating radiative pumping 
and location in the high radiation density environment of the nuclear regions. Abundance ratios OH + /H20 + of 1 - 2.5 are esti- 
mated in the nuclei of both sources. The inferred OH + column and abundance relative to H nuclei are (0.5 - 1) x 10' s cirr 2 and 
~ 2 x 10~ 8 , respectively. Additionally, in Arp 220, an extended low excitation component around the nuclear region is found to have 
OH + /H 2 + -5-10. H 3 + is detected in both sources with A'(H 3 + )~ (0.5 - 2) x 10 16 cm" 2 , and in Arp 220 the pure inversion, 
metastable lines indicate a high rotational temperature of ~ 500 K, indicative of formation pumping and/or hot gas. Simple chemical 
models favor an ionization sequence dominated by H + — > + — > OH + — > H 2 + — > H 3 + , and we also argue that the H + production is 
| most likely dominated by X-ray/cosmic ray ionization. The full set of observations and models leads us to propose that the molecular 

ions arise in a relatively low density (> 10 4 cm 4 ) interclump medium, in which case the ionization rate per H nucleus (including 
secondary ionizations) is f > 10~ 13 s , a lower limit that is several x 10 2 times the highest rate estimates for Galactic regions. In 
Arp 220, our lower limit for f is compatible with estimates for the cosmic ray energy density inferred previously from the supernova 
rate and synchrotron radio emission, and also with the expected ionization rate produced by X-rays. In NGC 4418, we argue that 
X-ray ionization due to an AGN is responsible for the molecular ion production. 

, Key words. Line: formation - Galaxies: ISM - Infrared: galaxies - Submillimeter: galaxies 

O 1 

in: 

^-H 1 ■ Introduction abundant species in the nuclei of some (ultra)luminous infrared 

i— H galaxies ((U)LIRGs). 

Observations of O-bearing molecular ions are powerful probes ~ xhe success ful launch and operation of the Herschel Space 

. of the oxygen chemistry and related physical processes in the observatory and its instruments have allowed the observation 

•• ; interstellar medium. The formation of OH\ H 2 + , and H 3 + of O-bearing molecular ions in both galactic and extragalac- 

. £ . can be i nitiated by cosmic ray and/or X-ray ionization of H and tic sources . i n the Milky Way, and foll owing the first d e tectio n 

X " H2 (e-g- iHerbst & Klempereij Mil |Maloney et al.j ll99fl), and of OH + towards Sgr B2 with APEX (IWvrowski et all [2010), 

^ ; hence their abundances are sensitive to the cosmic ray and X- the ground -state lines of OH + and H 2 + have been observed 

. CO . ray fluxes that permeate the observed environments. Significant in absorption toward a number of star-forming regions with 

amounts of molecular ions can also be formed in the surfaces high subm illimeter continua, mainly tracing intervining dif- 

of_rjhotodissociation regions (PDRs) (|Sternberg & Dalgarnoj f use /translucent c louds with low H? fraction (lOssenkonf et all 

[19951). The OH + /H 2 + and OH + /H 3 + abundance ratios are sen- 2 010; Ger inetdl 1201(1 iNeufeld et all l20ioh . The ground-state 

sitive to the molecular fraction and to the fractional ionization, lines of both species have also been detected in the Orion 

as the destruction paths for the three species are dominated by BN/ KL outflow, though still in absorption against the contin- 

dissociative recombination (H 3 Q + ) and reactions with H 2 (OH + uum dQupta et all EoiCh . and toward the high-mass star form- 

and H 2 + ; |Gerin et all |2010t |Neufeld et all |2Qlfl). In environ- ing regio n AFGL 2591, also ind icating low excitation for both 

ments with T < 250 K, these are the intermediate species for hydr i des (Bru derer et all iSal) . Most of the Galactic OH + and 

the main gas-phase production route of OH and H 2 0, which are Ul0 + detections essentially sample low-excitation gas and are 

thought to be generated by the cosmic ray field in our galaxy 

* Herschel is an ESA space observatory with science instruments pro- (OssenkopfluullpOirirGerin et all 12010: Neufeld et afl |2010|; 

vided by European-led Principal Investigator consortia and with impor- iHollenbach et all |2012|) . The more stable H 3 + was first de- 

tant participation from NASA. tected in the Sgr B2 molecular cloud complex near the Galactic 
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center (IWootten et al. [ 199l[). and fur t her m apped with APEX 
in the same region (Ivan der Tak et all l2006h through the non- 
metastable 3j - 2^ submillimeter line at 365 GHz. Recently, 
high-lying lines of H 3 Q + from me tastable levels have been de- 
tected in Sgr B2 bv lLis et al] d2012l) . who suggested that the ex- 
citation could be due to formation pumping in X-ray irradiated 
gas. 

The situation in luminous extragalactic sources is different in 
some aspects. While the ground-state tr ansition of ortho-H 2Q + 
has been detected in absorption in M82 dWeifi et all 12010). the 
ground-state OH + and HaO + lines have been detected in emis- 
sion in Mrk 231, probably indicating a combination of high 
excitation conditions and high ionization rates that are plausi- 
bly tracing X-ray Domin ated Regions (XDRs) due to the AGN 
( Ivan der Werf et aUl2010l) or the hi gh cosmic ray intensit ies sug- 
gested to be present in ULIRGs (Papadopoulosl I20T0I) . From 
ISO observations, hints of OH + absorption at ~ 153 /mi from 
the first excited rot ational level were previously re ported in the 
ULIRGs Arp 220 dGonzalez-Alfonso et all 12004 her eafter G- 
A04), and Mrk 231 dGonzalez-Alfonso et all 12008). but the 
absorption in Arp 220 was mostly ascribed to NH. Recently, 
Rangwal a et alj d201 ll) (hereafter Rl 1) have reported the detec- 
tion of the ground-state OH + and o-H20 + fine-structure com- 
ponents in absorption toward Arp 220, as well as two emission 
lines of H20 + with excited lower levels. On the other hand, the 
H 3 Q + 3t - 2z submillimeter line has also been detected in Arp 



220 dvan der Tak et all 2008, hereafter vdT08) and other nearby 



galaxies dAalto et all 1201 lb . 



We report in this work on the Herschel/PACS detection of 
excited OH + , H ? Q + , and H 3 Q + in both Arp 220 and NGC 4418 



Aid'. 

In a recent study dGonzalez-Alfonso et all 120121, hereafter Paper 
I), we have reported on the Herschel/PACS observations of these 
two sources in H 2 0, H' s O, OH, 18 OH, HCN, and NH 3 , together 
with a quantitative analysis in which a deconvolution of the con- 
tinuum into several components enabled us to derive the col- 
umn densities and abundances of these species associated with 
each of the components. A chemical dichotomy was found in 
Paper I: on the one hand, the nuclear regions of both galaxies 
have high H 2 abundances (~ 10~ 5 ), and high HCN and NH 3 
column densities, indicating a chemistry characterized by high 
Tgas, Tdust, and high density, where the release of molecules from 
grain mantles to the gas phase has occurred (through evapora- 
tion in "hot core" chemistry and sputtering in shocks). On the 
other hand, the high OH abundance found in both sources, with 
OH/H2O ~ 0.5, appears to indicate the occurrence of photo- 
processes due to X-rays and/or the effects of cosmic rays. This 
second aspect of the chemistry is further studied in this paper 
from the observed O-bearing molecular cations detected in the 
Herschel/PACS spectra of both sources. The paper is organized 
as follows: the observations are presented in Sect. |2j column 
densities and abundances are inferred in Sect. [3j chemical mod- 
els are presented and compared with observations in Sect. HJ the 
inferred important processes are discussed in Sect. [3] and the 
main conclusions are listed in Sect. [6] As in Paper I, we adopt 
distances to Arp 220 and NGC 4418 of 72 and 29 Mpc, and red- 
shifts of z = 0.0181 and 0.00705, respectively. 



2. Observations and results 

The full, high resolution, PACS spectra of NGC 4418 and Arp 
220, were observed as part of the guaranteed-time key program 



SHINING on July 27fh and February 27th (2010), respectively 
The data reduction was done using the standard PACS reduction 
and calibration pipeline (ipipe) included in HIPE 5.0 975. For 
the final calibration, the spectrum was normalized to the tele- 
scope flux and recalibrated with a reference telescope spectrum 
obtained from Neptun e observations. Furth e r deta il s of the data 
reduct ion are given in Gracia-Carpio et alJ (1201 lh . [Sturm et al. 
d2011l) . and Paper I. 

The PACS spectrum of Arp 220 in Paper I was corrected 
in both flux and wavelength calibration due to the effects of 
sub-spaxel pointing errors. Following newly made software cor- 
rections to the Herschel pointing star catalogue and the star- 
tracker x and y focal lengths, a second full spectrum of Arp 
220 was taken by the Observatory to assess the improvement to 
Herschel's pointing. Although the baseline of the new spectrum 
is of higher quality, general agreement with our corrected origi- 
nal spectrum was found in both flux and wavelength calibration. 
The data shown here are mostly an average of the two observa- 
tions, except for some ranges where the new data set yields a 
more reliable baseline. 

Spectroscopic parameters of OH + , H20 + , and H 3 + used for 
line identification and radiative transfer modeling were mostly 
taken from the spectral line catalogs of the CDMS dMiiller et all 
I2001L 120051) and JPL dPickett et all Il998l) . and are included in 
Table [1] together with the equivalent widths (W eq ) and fluxes of 
the detected lines. The wavelength uncertainties tabulated in the 
catalogs are significantly lower than the PACS spectral resolu- 
tion. Energy level diagrams for the three cations are shown in 
Fig. Q] where the solid red arrows indicate the detected transi- 
tions in at least one of the sources. 

2.1. OH + 

Owing to the electronic spin S = 1 , all rotational levels of OH + 
except the ground N - state are split into 3 fine-structure levels 
with total angular momentum / = N-l,N, N+ 1 (Fig.[T||3. Of the 
six allowed fine-structure transitions between excited rotational 
levels, the three with AJ = AN are the strongest. Line spectro- 
scopic par ameters of OH + wer e taken from the CDMS, and were 
derived bv lMiiller et al.ld2005l). In that wor k, theN - 1-0 transi- 
tion frequencies from Bekoo v et al.l(ll985l) were fit together with 
additional infrared data. A dipole moment o f 2.256 D was taken 
from an ab initio study (Werner et al., 1983). 

Figure|2]shows the observed spectra around the OH + lines in 
NGC 4418 (upper histograms) and Arp 220 (lower histograms). 
The spectra around the positions of the 2j <— l /< lines, with 
Slower ~ 50 K, are displayed in panels a-c. The strongest 2 3 <— I2 
at w 153 fim is partially blended with NH, but the 22 * — 1 1 is 
free from contamination and clearly detected in both sources. 
In Arp 220, the 153 fim absorption appears to be dominated 
by OH + rather than by NH, in contr ast with the initial esti- 
mate by iGonzalez- Alfonso et alJ (T2004) based on ISO data. The 
2 1 <— lo at » 148.7 fim is also well detected, and the intrinsi- 
cally weaker 2i<— 1] (158.4 yum) and 22 <— I2 (147.8 //m) lines 
are (marginally) detected in NGC 4418 but still very strong in 
Arp 220. There is in addition an apparent red-shifted wing in 
the NH 2 3 <— I2 spectrum of Arp 220, detected in both data sets, 



1 The observation identification numbers (OBSIDs) are 
1342202107 - 1342202116 for NGC 4418, and 134191304 - 
1342191313 and 1342238928 - 1342238937 for Arp 220. 

2 The hyperfine structure of both OH + and H20 + is ignored in this 
work as the line splitting is much lower than both the PACS spectral 
resolution and the linewidths. 
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Fig. 1. Energy level diagrams of OH + , H20 + , and H30 + , showing with arrows the transitions that lie within the PACS range; colored 
numbers indicate round-off wavelengths in yum. The fine-structure splitting of the OH + and H20 + levels is out of scale. Solid-red 
arrows mark the lines detected in NGC 4418 and/or Arp 220; dashed-red arrows indicate marginal detection in any of the sources; 
green arrows indicate blended lines, and black arrows show undetected lines. 



which coincides with the position of the intrinsically weak OH + 
2i <— 12 transition. As shown below (Sect. 12.31 ). however, this 
wing feature is most probably due to the very excited 12^ <— 
12+ 2 line of H 3 + . 

Most 3 j <— Iji lines at ~ 100 //m lie within the gap between 
the green (A < 100 //m) and red (A > 100 //m) bands of PACS 
and are not observable, except for the 32 <— 22 transition at 103.9 
fj.m shown in Fig. [2]J. The line in NGC 4418 is close to the edge 
of the red band where the spectral noise is high, but in Arp 220 a 
clear spectral feature is detected at the expected wavelength. The 
line, however, is blended with the H2O 615606 transition (Paper 
I), which may account for some of the observed absorption. 

The three strongest 4y <— 3 j< lines at 76 - 76.5 yum are clearly 
detected in both sources (Fig. |2^), indicating a high-excitation 
component in the transient OH + . It is worth noting that while the 

2 j <— I j> lines are stronger in Arp 220, the more excited 4j <— 

3 j> lines are slightly deeper in NGC 4418, resembling the behav- 
ior found for other species like H2O, OH, HCN, and NH3 (Paper 



I). The 4j <— 3j< lines still have slightly higher equivalent widths 
in Arp 220 (Table []]) owing to the broader line profiles in this 
source. The line opacities of the fine-structure components are 
r oc A 3 g u A u [ if the lower levels are populated according to their 
degeneracies, and then T4 5 _3 4 : T4 4 _3 3 : T4 3 _3 2 = 1 : 0.77 : 0.58. 
These are expected to be the ratios of the line equivalent widths 
(Weq) if the three components are optically thin. In Arp 220, the 
observed W eq ratios are 1 : (0.69 ±0.11) : (0.53±0.10)(TableQ), 
consistent with optically thin absorption. In NGC 4418 they are 
1 : (1.2 ± 0.6) : (1.0 ± 0.5), compatible within the uncertainties 
with both optically thin and saturated absorption. 

In NGC 4418, some marginal features around 61 fim roughly 
coincide with the positions of the OH + 5/ <— 4j> lines (Fig. [2}). 
However, the feature associated with the intrinsically strongest 
56 <— 45 component is relatively weak and shifted in wavelength 
by 0.015 /im, shedding doubt on this tantalizing identification. In 
Arp 220, this part of the spectrum was affected by instrumental 
problems and is not shown. 



3 



Gonzalez-Alfonso et al.: Excited OH + , H 2 + , and H 3 Q + in NGC 4418 and Arp 220 



o 

0) 

a 0.8 

m 

CD 
N 



0.6 



a) 2 g - 


i i i i i i 

1 O 

-lj 2 3 - 


1 Q 


1 


b)2 s 


i i i i i i i i i 

1 




i 

cj 2 r 


1 l"Ln 














r- 











/A 

7 ^ 

jj/jv NH | 


H 3 + - 
. 




^ \ 




T- 


[ ^ 
i 


f - 

x2 - 



u 
o 

a 

a 
'■+-> 

a 
o 
u 



152 



152.5 



153 



153.5 



148 



148.5 



158.4 158.8 



0.95 



;; e ) 4 3 -3 2 4 4 -3 3 4 5 -3 4 ;!f) 5 4 -4 3 5 5 -4 4 5 6 -4 5 




103.5 104 



76.2 



A 



rest 



76.4 
(yum) 



76.6 



61.2 



11.4 



Fig. 2. Continuum-normalized spectra around the OH + lines in NGC 4418 (upper histograms) and Arp 220 (lower histograms). The 
dotted vertical lines indicate the rest wavelengths of the transitions for the nominal redshifts of z = 0.00705 (NGC 4418) and 0.0181 
(Arp 220). The spectra in panels b and c have been scaled by a factor 2. (Potential) contamination by other species (NH in panel a 
and H2O in d) is also indicated. The marginal 5j -4j> lines in f correspond to NGC 4418; this part of the spectrum is not available 
in Arp 220. Blue lines show model fits for OH + in the nuclear components of both NGC 4418 and Arp 220, the green lines show 
the model fit for OH + in the extended components of both sources, and the light-blue lines in a are model fits for NH; red is total. 



2.2. H 2 + 

The asymmetric H20 + also has fine-structure splitting of the ro- 
tational levels producing doublets (S = 1/2) and yielding a very 
rich spectrum in the far-infrared domain (Fig. Q]). Predictions 
of the rotational spectrum were taken from the CDMS. The 
C DMS data w e re bas ed on extrapolated zero-field frequencies 
in iMurtz et all d 19981). the derived jV =1-0 transition fre- 
quencies based on lStrahan et al. (1986), and additional rovibra- 
tional (infrared) transitions and ground state data from elec- 
tronic spectra. HIFI observations of the = lio - lot transi- 
tion frequenc i es abo ve 600 GHz toward Sgr B2(M) reported by 
ISchilke et al.l (1201 Oh are compatible with the transition frequen- 
cies in the CDMS, although these lines lie outside our accessible 
frequency range. A gr ound state dipole moment of 2.398 D was 
calculated ab initio bv lWeis et al.l dl989l) . 

As shown in Fig. [T] most low-lying lines of H20 + are 
blended with lines of other species, but some of them lie in 
free windows (Fig. |3). Very prominent is the ground-state p- 
H 2 + 212101 transition at 183 - 184 /mi in Arp 220, with the 
three fine-structure components clearly identified in the spec- 
trum. This transition is not detected in NGC 4418. In Arp 220, 
the observed line strengths strongly decrease with increasing 
£)ower, indicating that the ground-state p-H20 + lines arise in 
low excitation gas. The opacity ratios of these components are 



T5/2-3/2 : T3/2-1/2 : T3/2-3/2 = 1 : 0.55 : 0.12, so the observed 
W eq ratios of 1 : (0.65 ± 0.07) : (0.30 + 0.05) (Table [Q suggest 
line saturation effects in at least the strongest 5/2-3/2 com- 
ponent. For cr v = 60 kms -1 , this component attains an optical 
depth of unity with a column of iV(p - H20 + ) ~ 5 x 10 14 cm -2 . 
The line absorbs only w 15% of the continuum, likely indicating 
some spectral dilution and a low covering factor of the contin- 
uum. 

The two prominent features at 104.7 - 105 yum in both 
NGC 4418 and Arp 220 are most probably due to NH 2 . The fea- 
ture at 143.8 fim may be contaminated by C3. In Arp 220, the ex- 
cited 221 1 10, 3"13202, 32221 1, and 404313 lines of H 2 0\ with 
Slower up to 190 K, are clearly detected, but in NGC 4418 only 
the high-lying 322211 transition show measurable absorption. 
Detection of the 404313 transition in NGC 4418 is question- 
able, as the spectral features associated with the 9/2 - 7/2 and 
7/2 - 5/2 components are either too narrow or shifted (Fig. [3^). 

2.3. H 3 + 

H30 + and NH3 are isoelectronic and spectroscopically similar: 
both are oblate symmetric top rotors, have a double-minimum 
potential such that inversion takes place between the two equiva- 
lent minima, and have ortho (K - 3n) and para (K + 3n) modifi- 
cations. Fig.Q]shows the energy level diagram for H30 + . Owing 
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Fig. 3. Continuum-normalized spectra around the H20 + lines in NGC 4418 (upper histograms) and Arp 220 (lower histograms). 
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to the AK = selection rule for all radiative transitions, the 
J - K levels in the lower tunneling state cannot be radiatively 
pumped and hence are metastable. They therefore will be "ther- 
mally" populated at a rotational temperature, T mt , which reflects 
either r gas or, if the formation-destruction rates are faster than 
the rate for col lisional relaxation, the formation pro cess char- 
acteristics (e.g. iBruderer et ail l2010bt iLis et all 120121 see also 
AppendixIBl. 

Spect roscopic da t a for H3Q + were taken from the JPL cat- 
alogue (Pickett et al., 1998). The entry was based on the recent 
analysis bv lYu et all d2009l) . who also presented some improved 
transition frequencies. The majority o f the pure tunneling t ran- 
sition frequencies were reported by Verh oeve et al.l dl989h . A 
gr ound state dipole momen t of 1 .44 D was calculated ab initio 
bv iBotschwina et aD (Il985|). The antisy mmetric^ state is higher 
in energy by 1659 GHz (lYu et al.U 2009> because of the low bar- 
rier to linearity, and thus the pure inversion transitions of HaO + 
lie in the far-IR domain (150- 185 //mfj 

The solid-red arrows in Fig. Q] indicate the HsO + lines de- 
tected in Arp 220, and Fig.|4]presents the observed H30 + spec- 
tra of NGC 4418 and Arp 220. In Arp 220, spectral features at 



3 The symmetric and antisymmetric states are here denoted with + 
and -, respectively. 

4 For NH 3 , the antisymmetric state is higher in energy by 23.7 GHz 
and the pure-inversion transitions lie at centimeter wavelengths. 



the wavelengths of the H30 + pure-inversion, metastable (here- 
after PIMS) lines (panels a-i) are detected up to a lower level 
energy of * 1400 K above the ground state (12 + 2 ), while none 
of these lines is detected in NGC 4418. The rotation-inversion 
(hereafter PJ) lines lie at shorter wavelengths (A < 100 yum) and 
are shown in Fig.|4j-p- Four of them arise from non-metastable 
levels (Fig. [TJ, three of them from metastable levels (hereafter 
RIMS lines), and are detected within the lowest K < 4 ladders 
up to an energy of ~ 200 K in both galaxies. 

The l~ - 1|, 2" - 2 + , and 11" - 11+ PIMS lines are con- 
taminated by strong lines of CH, H^O, and [C n], and are not 
shown. All features in Arp 220 are detected in both datasets. 
Close lines due to other species, as well as potential contamina- 
tions (panels f and i), are also indicated in Fig. [4] Based on the 
weakness of the HCN 19-18 line (Paper I), the HCN 20-19 
line at 169.39 fim is not expected to strongly contaminate the 
H30 + 8 8 <— 8 + line, though some contribution from HCN is not 
ruled out (see below). Likewise, the nearby NH2 330321 4 -3 
line is intrinsically weak, and NH2 modeling that accounts for 
other detected NH2 lines does not produce significant absorption 
at 169.38 fim. A broad feature is detected at 165.7 fim (panel 
g), which cannot be explained by only the NH3 3 2 <— 2 + and 
3^ <— 2| lines (Paper I); thus strong absorption by ortho-H30 + 
9 9 <— 9 + is inferred. The relatively weak H30 + 10[ <— 10 + line 
is affected by an uncertain baseline, though the line is detected 
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in both data sets. On the other hand, a line wing feature close 
to the NH 23 <— I2 153.35 fim line is detected, coinciding with 
both the ortho-H 3 + 12^j e- 12[ 2 and the OH + 2 t <- 1 2 lines 
(153.47 fim, see also Fig.[2^). As mentioned earlier, the latter is 
also intrinsically weak, and modeling of the other OH + 2j <— 1 /< 
lines (Sect. 13.4b yields negligible absorption in that component. 
Upon inspection of the CDMS and JPL catalogs, we have not 
found any other reliable candidate for the NH wing, which we 
tentatively ascribe to the very high-lying H30 + line. 

No potential contamination has been found for the RI lines. 
Contrary to the PIMS lines, some of these RI lines are detected 
in NGC 4418. The non-metastable levels are populated through 
absorption of far-IR photons, and trace the local far-IR radiation 
density. 

The PIMS and RIMS lines are observed in absorption against 
the strong far-IR continuum source of Arp 220, and their relative 



strengths probe T vot . The population diagram of H30 + is shown 
in Fig. [5] where the level columns have been calculated in the op- 
tically thin limit and by assuming that the H30 + covers the whole 
continuum from Arp 220. Since the degeneracy of the ortho lev- 
els is twice that of the para levels, we expect a nearly continuous 
distribution for an ortho-to-para ratio of 1, as observed in the 
high-lying levels (crosses in Fig. [5]). The distribution is roughly 
described with two values of r rot , ~ 180 K up to the 5+ level, 
and ~ 550 K for higher-lying l evels. This is strikingly similar to 
that found bv iLis et aD (1201 2h towards Sgr B2. The column of 
the 8g level (E/ as 660 K) appears to be higher than expected, 
suggesting some contamination by HCN at 169.4 fim (Fig.|4f). 
On the other hand, there are two metastable levels, the 3J and 
the 4J, for which both PIMS and RIMS lines are available. In 
both levels, the columns derived from the RIMS lines are higher 
by a factor of a; 1.8 than those derived from the PIMS lines. 
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Two main reasons can account for this: first, the level columns 
are calculated, as mentioned above, by assuming that the H30 + 
covers the whole far-IR continuum, which is relatively cold (i.e. 
S i8o/^60 = 0.43). If the H30 + is absorbing a warmer continuum 
that is diluted within the observed SED (i.e. filled in by a cooler 
continuum source), the intrinsic equivalent width and column of 
the PIMS lines will be increased relative to those of the RIMS 
line, and both columns will be in better agreement. Therefore, 
the discrepancy probably indicates that the highly excited H30 + 
covers only a fraction of the observed continuum, and that the 
level columns are lower limits. Second, the columns in Fig. [5] 
are not corrected for stimulated emission, which will mainly in- 
crease the columns of the PIMS lines at long wavelengths. 

2.4. Velocity shifts in Arp 220 
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Fig. 5. Population diagram of H30 + (metastable levels) in 
Arp 220. Red and blue symbols correspond to ortho and para 
levels, respectively. Columns derived from pure-inversion and 
rotation-inversion lines are shown with triangles and squares, 
respectively. The crosses indicate the expected population dis- 
tribution for r iot = 550 K and a fixed ortho-to-para ratio of 1 . 



Having confirmed the velocity calibration (or "corrections") 
applied to the Arp 220 spectrum (Paper I) with the second Arp 
220 dataset, we now compare the line profiles and velocity shifts 
presented in Paper I to the velocity profiles and centroids of 
the molecular ions. The average profiles of several H2O, OH + , 
H 2 0\ H 3 + , HCN and NH 3 lines detected in Arp 220 are com- 
pared in Fig. [6] In all panels, HE and LE mean high-excitation 
and low-excitation, respectively, qualitatively indicating whether 
the selected lines are high-lying or low-lying lines. Thus H2O 
HE is the average profile of two representative high-excitation 
lines of H 2 0, the 616505 and 322211 lines (Paper I); HgO LE 
is the average profile of the two low-lying 212101 and 221110 
lines, and similarly for other species/excitation. The H30 + PIMS 
spectrum is the average of the unblended high-lying, pure- 
inversion metastable lines of H30 + , and the NH3 RIMS one is 
the average of the relatively uncontaminated rotation-inversion 
metastable 7£ <— 6 6 and 6 5 <— 5g lines. The averages were gen- 
erating by resampling the individual spectra to a common veloc- 
ity resolution, and adopting equal weights for all them. Figure|7] 
displays the velocity centroids of these generated line profiles. 



Figures [6] and [7] indicate that there are systematic velocity 
shifts associated with line excitation, irrespective of the con- 
sidered species. The averaged H2O HE profile, peaking at -20 
km s" 1 , has a shape very similar to the averaged OH + HE (panel 
a), H 2 + HE (b), and NH 3 HE (d) profiles, though the latter is 
somewhat narrower (and relatively uncertain due to line blend- 
ing). The HCN 18 <— 17 line (panel a) could be even more 
blueshifted, though the signal-to-noise ratio is relatively low. On 
the other hand, there is a clear velocity shift in the peak absorp- 
tion of the H2O HE and LE averages, with the latter peaking at 
a velocity of +50 kms -1 . The LE averages of OH + and H20 + 
also tend to peak at more positive velocities, though the OH + LE 
lines are not ground-state lines CEi 0W er ~ 50 K) and thus may 
represent an intermediate situation. 

The velocity centroids of the H30 + PIMS lines, however, 
which are characterized by high Ei OW er, peak at the velocity of 
the LE lines and thus present an exception to the finding dis- 
cussed above. A velocity shift of as 70 km s is shown in Fig.[6j; 
between the average H2O HE and H30 + PIMS spectra, with the 
latter peaking at the velocity of the LE lines of H2O and H20 + . 
This shift is much lower than the line widths of ~ 300 kms -1 , 
but larger than our spectral calibration accuracy. The average ve- 
locity of all H2O lines (excluding the lowest-lying LE lines) is 
-20 + 35 kms -1 , whereas that of the metastable H30 + lines is 
50 + 25 kms -1 . The latter corresponds to vlsr = 5480 kms -1 , 
in rough agreement with the velocity centroid of the 3 J — 2^ line 
(vdT08) an d close to the central velo city of CO 2-1 in the west- 
ern nucleus dDownes & Eckartll2007l) . Excluding the low signal- 
to-noise H30 + 4^ <— 3j line, the other three non-metastable 
lines have an average velocity of 20 + 15 kms -1 . The low-lying 
4j <— 3j line (Fig. |4]s;) is even more redshifted (« 100 kms -1 ), 
with the two data sets giving very similar line shapes. With the 
same lower level, the pure-inversion 3^ <— 3^ line (Fig.|4^) also 
shows hints of redshifted absorption at > 100 km s -1 . 

At the very least, the velocity shift between the peak absorp- 
tion of the H30 + PIMS and the HE lines of other species in- 
dicates a spatial shift between the regions responsible for the 
strongest absorption in both sets of lines. On the one hand, this 
may not be surprising from excitation arguments. The H2O HE 
levels are excited through radiative pumping by dust (Paper I), 
and thus the corresponding lines trace the region of highest far- 
IR radiation density (and thus of highest Tjust)- On the other 
hand, the H30 + metastable levels cannot be radiatively pumped, 
but are either excited through collisions indicating high r gas , or 
reflect H30 + formation in high-lying levels; thus the physical 
conditions that favor strong absorption in each set of lines are 
different and do not necessarily coincide spatially. On the other 
hand, the NH 3 RIMS lines, with £ lower = 300 - 400 K, are colli- 
sionally pumped but better match the H2O HE profile (Fig. |6}l), 
suggesting either chemical differences between the two veloc- 
ity components, or the dominance of chemical pumping on the 
formation of the PIMS lines of H 3 + . 

Given the velocity coincidence between the peak absorption 
of the H30 + PIMS lines and the H2O LE lines, the question is 
whether the H30 + PIMS lines arise from foreground gas in the 
direction of, but detached from the nuclear region, or the lines 
arise on the contrary from gas physically associated with the nu- 
clei. The PACS spectrum of Arp 220 shows RI lines of H30 + 
arising from non-metastable levels (Fig. |4jn,l,n,p) that are ra- 
diatively pumped through dust emission. Of the above lines, the 
H30 + 4 <— 3q at 70.8 /im is most sensitive to rd ust , requiring 
relatively high far-IR radiation densities. Figure [6}i indicates a 
good match between the shape of the H30 + 70.8 pm line and of 
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Fig. 6. Comparison between averaged line shapes of several lines in Arp 220. HE, LE, PIMS and RIMS mean high-excitation, 
low-excitation, pure-inversion metastable, and rotation-inversion metastable, respectively. H2O HE is the average of the 616505 and 
322211 lines of H 2 0(Paper I); HCN HE is the 18 - 17 line of HCN; OH + HE is the average of the 4 3 - 3 2 and 4 5 - 3 4 lines of 
OH + ; H 2 + HE is the average of the 322211 7/2 - 5/2, 313202 5/2 - 3/2, and 404313 7/2-5/2 lines of H z O + ; H 3 + PIMS is 
the average of the pure-inversion, metastable 5 5 - 5t, 6 6 - 6£, 7^ - lij, and 8g - 8^ lines of H30 + ; NH3 RIMS is the average of 
the rotation-inversion metastable 7£ - 6 6 and 6<r - 5 J lines of NH3 ; H 2 LE is the average of the 212101 and 22 1 1 1 lines of H 2 0; 

OH + LE is the average of the 2 2 - 1 2 , 2 X - 1 , 2 2 — 1 1 , and 2 X - U lines of OH + ; H 2 + LE is the average of the 212101 5/2 - 3/2 
and 3/2 - 1/2 lines of H 2 + ; H 3 + 70.8 is the 4 Q - 3q line of H 3 + at 70.8 fim. The dotted and dashed vertical lines indicate the 
velocity centroids of the averaged HE and LE spectra, respectively. 



the PIMS average, and thus it is likely that both sets of lines arise 
from essentially the same region. We thus expect that the H 3 + 
PIMS lines are of primarily nuclear origin though, as argued in 
Sect. 13.4.31 probably more extended than the region where the 
NH 3 lines are formed. 

The simplest scenario that can account for the observed ve- 
locity pattern in Arp 220 consists of the H 3 + PIMS lines and 
the LE lines of other species peaking at the "rest" velocity of 
the western nucleus, and thus primarily tracing rotating gas to- 
gether with more foreground gas, while the HE lines of H 2 0, 
OH, NH 3 , and also OH + and H 2 + , mainly tracing an outflow- 
ing component with a line of sight velocity of ~ 75 km s _1 , sim- 
ilar to that observed in HCO + (Sakamoto et al. , 2009fl Besides 
grain mantle evaporation and neutral-neutral reactions at high 
r gas (Paper I), sputtering in shocks could also remove the man- 
tles and contribute to enhance the gas-phase abundances of H 2 
and NH 3 . Alternatively, both components may be tracing differ- 
ent sections of the nuclear rotating disk(s), which shows steep 
velocity gradients as seen in the (sub)millimet er lines of CO 
(Downes & Eckartl 120071; ISakamoto et aUl2008l) . 
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Fig. 7. Velocity centroids of the averaged line shapes in Arp 220 
displayed in Fig. [6] derived from Gaussian fits. Velocities are 
calculated relative to z = 0.0181. Blue and red symbols corre- 
spond to HE and LE averages, respectively; for H 3 + the green 
symbols indicate the PIMS and RI averages. Errorbars are lcr 
uncertainties. 



5 In Fig. 12 of Paper I, the HCO + lines should be shifted 100 km s _I 
to the red, as the original spectrum 1 Saka moto et all 12009) uses the ra- 
dio convention for the velocities; thus the HCO + absorption toward the 
western nucleus better matches the blueshifted absorption in the H 2 
HE lines. 



3. Environments, column densities and abundances 

3.1. The excitation of the molecular ions 

Derivation of the column densities implied by the observations 
require calculations of excitation appropiate for the environ- 
ments where OH + , H 2 + , and H 3 + reside. The rates for colli- 
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sional excitation of these molecular ions are not known, but both 
Arp 220 and NGC 4418 have luminous-compact nuclear far-IR 
continuum sources that excite H2O, OH, and other species to 
level energies of > 500 K (Paper I). Therefore, a natural mecha- 
nism for exciting OH + , H20 + , and the non-metastable levels of 
H30 + along a given /f-ladder is via radiatively allowed transi- 
tions through absorption of dust-emitted photons. 

Collisional excitation of the high-lying lines of OH + and 
H20 + is expected to be less efficient. The A„/-Einstein coeffi- 
cients are high (TableQ}, and thus the critical densities (n cl ) are 
likely to be high as well. For the N - I - transition of OH + , 
with A„i =s 0.02 s _1 and assuming a coefficient for collisional 



for formation pumping to dominate the observed excitation over 
collisions can be written a^| 



de-excitation with H and H2 of 10 cm s 



2 x 10 7 



cm , which is much higher than the expected densities (see 



Sect 14.5 1 below), and also higher than the nuclear densities in- 
ferred from other tracers (Paper I). As a result, the expected ex- 
citation temperature of the N = 1 - transition of OH + from 
collisional excitation alone would be < 20 K, in contrast with 
the typical values of 35 - 70 K obtained from radiative exci- 
tation in the models shown below. On the other hand, excita- 
tion through collisions with electrons has rate coefficients that 

10~ 6 cm 3 s" 1 



migh t be much highe r than with H a nd H? , 
dNeufeld & Dalgarno .ll989t iLim et all 1 1999h . and has been re 



cently invoked bv Ivan der Tak et al.l (120121) to explain the emis- 
sion in the HF 1-0 line toward the Orion Bar. However, the 
expected electron density is several x 10 3 lower than the density 
of H (Fig. IC. Il l, so that its effect on the excitation is not expected 
to be much higher than that of H. Finally, the detected lines are 
seen in absorption, and radiative pumping models generate the 
observed absorption features naturally, as we show below. 

Another potential excitation mechanism is c hemical forma- 
tion in excited levels (IStauber & Brudereiil2009l) . i.e. formation 
pumping. In environments with weak dust emission, a down- 
ward cascade after m olecular formation in excited levels (e.g. 
Brud erer et all 1201 Obi) is expected to generate emission in non- 
metastable lines (i.e. chemoluminiscence). But again, in the spe- 
cific case of NGC 4418 and Arp 220, all detected lines from 
molecular ions are seen in absorption, qualitatively matching 
the far-IR pumping scheme. Furthermore, ignoring excitation 
by dust emission, the equilibrium level populations of OH + and 
H2Q + resulting from formation pumping (see Appendix B in 
iBruderer et all 1201 Oal) yield very low excitation even for very 
high formation rates (see Appendix lAb. This is due to the very 
high A-Einstein coefficients (Table [1]) as compared with the ex- 
pected formation/destruction rates (Sect. [4]). Thus we do not ex- 
pect that formation pumping has a measurable effect on the ex- 
citation of non-metastable transitions in these sources. 

The case of the metastable lines of H30 + is different. The 
K > 1 ladders cannot be populated from K - 0,1 through ra- 
diative pumping, but only through collisions or via molecular 
formation in high-lying levels. The latter process can be envi- 
sioned as follows: molecular formation in an excited (K, J) level 
will be followed by quick cascade down to the metastable J — K 
one, where no more paths for radiative decay are available. Then, 
the molecule will "freeze-out" in that metastable level (ignoring 
here radiative excitation for simplicity) until destruction through 
recombination or until a collisional thermalization event takes 
place. If the destruction rate dominates over collisional relax- 
ation, the observed population distribution in the metastable lev- 
els (i.e., riot) will be governe d by the formation process rather 
than by r gas dLis et all l2012h . Writing the rate of H 3 + forma- 
tion per unit volume as Fh 3 o+ = eH 3 o+«H£, where eH 3 o+ is the 
efficiency with w hich ionizations are t ransfered to the produc- 
tion of H 3 + (see lNeufeldet all 12010. for OH + ), the condition 



ir ,(T ) < 8 x 10- 11 x ^ x 

/CcolUgas)<8XlO X QA X 1() _ 17cm 3 



5x10-8 «N 



where k co i (in cm 3 s _1 ) is the rate for collisional relaxation from 
a metastable level at T gas , £ is the ionization rate per H nu- 
cleus, x denotes the abundance relative to H nuclei, and the ref- 
erence value of 6h,o+ = 0.4 is based on the chemical models 
of Sect. |4] Collisio nal rates for H 3 Q + e xcitation with H2 have 
been calculated by lOffer & van Hemertl ([19920 for a tempera- 
ture of 100 K and for the lowest H 3 + levels in the K = 0, 1 , 2, 3 
ladders. Taking a s a guide k c0 \ ~ 2 x 10~ 10 cm 3 s at 100 K 
from the 33 level (lOffer & van Hemertl fl9921. formation pump- 
ing may have important effects at high £/nn, but collisional rates 
for higher levels and with H2 and H are required to check this 
point. Calculations for excitati on of H 3 + through electron im- 
pact (Faure & Tenn vsonl 120031) indicate that the AK > transi- 
tions are negligible, so collisions with electrons are probably not 
important in populating the K > 1 metastable levels. 

In ou r calculations for H 3 Q + bel ow, we use the de-excitation 
rates bv lOffer & van Hemertl (Tl9921. and correct the excitation 
rates according to the value of r gas in our models and the re- 
quirements of detailed balance. Since no collisional rates that 
connect the low- and high-TT ladders are available, we simply 
force the metastable high-lying Jj* levels to be populated ac- 
cording to the adopted T mt , which represents r gas , the signature 
of the formation process, or a combination of both. This will be 
further discussed in Sect. [5] Additionally, we attempt to estab- 
lish in Appendix [B] under which conditions the formation pro- 
cess can generate the observed high excitation in environments 
with moderate T gas . Excitation of non-metastable levels is dom- 
inated by the far-IR field, and significant discrepancies between 
observations and model results would indicate the relevance of 
alternative excitation mechanism. 



3.2. Source components 

As shown in Paper I, the far-IR spectra of NGC 4418 and Arp 
220 cannot be described by a single set of ISM parameters, but 
different lines have different excitation requirements and are thus 
formed in different regions of the galaxies associated with differ- 
ent far-IR continuum components that we summarize here (see 
also Figures 1 & 16and Tables 1 & 2 in Paper I). For NGC 4418, 
at far-IR wavelengths we require (z) a nuclear core (C core ) with 
Tdust ~ 130 - 150 K, that provides absorption in the high- 
lying lines of H2O, OH, HCN, and NH 3 ; (//) a warm component 
(C wann ) with 7d us t ~ 110 - 85 K, that provides absorption in 
moderately excited lines of H2O and OH as well as a significant 
fraction of the far-IR continuum emission; and (Hi) an extended, 
possibly infalling component (C ex tended) with 7d us t ~ 90 - 40 K, 
which accounts for low-lying redshifted lines of OH and O . For 
Arp 220, we require (z) the western and eastern nuclear compo- 
nents (C west and C eas t), where the high-lying molecular lines are 
formed; (z'z) the extended component (C ex tendedX which provides 
a significant fraction of the far-IR continuum emission (though 
mostly reemission from the nuclei) and accounts for the ab- 
sorption in moderate and low-excitation lines of H2O and OH; 



6 In steady state, this condition can also be expressed as k m \ < 
lQ- i0 x (7^/lOOr 05 x 0r(e-)/1.4 x 1(T 4 ) cm 3 s _1 , suggesting possi- 
ble significant formation pumping effects in high-Ay translucent clouds 
(see Appendix iBt. 

7 Available at http://www.strw.leidenuniv.nl/~moldata/ 
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and (Hi) an additional low-excitation, absorbing component or 
"halo" (ChaioX with no associated continuum that is required to 
fit the absorption in the ground-state lines of H2O, OH, and O . 

Calculations of excitation, line fluxes, and line profiles were 
first performed for the nuclear components in both sources to fit 
the high-lying lines, and then were done for the extended compo- 
nents to complete fitting of the absorption in the low-lying lines. 
Table |2]gives the inferred column densities, based on a "screen" 
approach for the nuclear components (i.e. the molecules are lo- 
cated in front of the nucleus) due to the high continuum opacities 
of these components in the far-IR; the "mixed" approach (i.e. 
the molecules and dust are mixed) would give columns h ighly 
exceeding the listed values (iGonzalez-Alfonso et all 12008). For 
the Cextended, the "mixed" approach is used owing to the much 
more moderate continuum opacities (Paper I). 

3.3. NGC 4418 

3.3.1. OH + 

In NGC 4418 we first used the parameters of the C core and C warm 
components as given in Paper I to match the absorption in the 
high-lying 4j <— 3j* lines of OH + , by assuming a covering fac- 
tor of 1 . The OH + column densities were similar for both compo- 
nents, JV(OH + ) = (0.6 - 1) x 10 16 cirT 2 , respectively. The model 
fits were also very similar for both components, so that we are 
not able to determine in which nuclear component the OH + lines 
are (primarily) formed. By using the same normalization column 
as in Paper I, iV(H) = 4 x 10 23 cm~ 2 , the OH + abundance in the 
nuclear region of NGC 4418 is *(OH + ) = (1.5 - 2.5) x 10~ 8 . 
Contrary to observed, the 4^ <— 34 component in Fig. [2^ is pre- 
dicted to be the strongest, which could indicate significant de- 
partures of the fine-structure distribution of populations relative 
to the predicted one. The effective spatial extent (diameter) pre- 
dicted by the C core and C war m models is * 20-30 pc, respectively, 
andr dust = 150- 110K. 

Figure [2] shows (blue curves) the model fit for the nuclear 
component (similar for C core and C warm ), together with a model 
for NH in panel a (light-blue curve). The nuclear model that 
fits the 4 j <— 3y< lines underpredicts the 2% <— lj line and the 
153.0 //m feature (composed of the OH + 23 <— I2 and the NH 
22 « — 1 1 lines, see Fig. |2^). The OH + 2%—\\ line has a red- 
shifted excess, and the 153.0 //m feature cannot be fully repro- 
duced at just the central wavelength of the feature, which is also 
redshifted relative to the OH + line. This strongly suggests that 
the absorption that cannot be fit with the nuclear model is due 
to C ex tended, a more extended, inflowing component proposed in 
Paper I in order to explain the redshifted absorption observed 
in the ground-state lines of OH and [O 1]. Fitting the remaining 
OH + absorption with the C exte nded model in Paper I (green curves 
in Figure 2a-upper), the fit shown with the red curve in Figure 
2a is obtained. We infer jV(OH + ) ~ 3 x 10 15 cm" 2 , and roughly 
X(OH + ) ~ 2 x 10~ 8 . 

In summary, the OH + lines are essentially formed in the nu- 
clear C warm + C cole components, with N(OH + ) « (6 - 10) x 10 15 
cm 4 . Some contribution to the 2i <— \\ and 23 <— I2 lines due to 
inflowing gas in C ex tended is also inferred, with A^(OH + ) ~ 3xl0 15 

-2 

cm . 

3.3.2. H 2 + 

H z O + is weak in NGC 4418 (Sect.O, though the 5/2 <- 3/2 
component of the 322211 transition at 89 //m, and more tenta- 
tively the 7/2 <— 5/2 one, are detected in the source (Fig.[3]I). 



We have carried out models for the C core and C warm components 
that fit the above lines. Both models give similar columns of 
Af(H 2 + ) ~ 3 - 4 x 10 15 cirT 2 , but fail to account for the some- 
what marginal 404313 components (Fig. [3^). The C core model 
prediction better matches these features but, as shown with blue 
curves in Fig. [3p, may (marginally) overestimate the 313202 
components. The H20 + 331220 7/2 <— 5/2 at 65.6 fim (panel 
f) is blended with NH2 and 18 OH lines and affected by an uncer- 
tain baseline (Paper I). We conclude that the OH + /H20 + column 
density ratio in the nuclear region of NGC 4418 is in the range 
1.5 — 2.5. The ground-state H20 + lines are not detected, indicat- 
ing that H20 + has a low abundance in the C ex tended component of 
NGC 4418. 

3.3.3. H 3 + 

Only the RI lines of H 3 + were detected in NGC 4418. They 
were fitted with the C core and C waim components, with an as- 
sumed ortho-to-para ratio of 1. Both components give similar 
column densities of jV(H3Q + ) *_5 — 8 x 10 15 cm 2 , similar 
to that inferred by Aalt o et alj d201 lb from observations of the 
3^ <— 2 2 line. The C core model is shown with blue curves in 
Fig. [4] The main drawback of both models is that the absorp- 
tion in the 3^ <— 2| transition (panel p) is significantly overpre- 
dicted; the line is not detected, possibly due to cancellation be- 
tween absorption toward the continuum source and reemission 
from a more extended region. The OH + /H30 + column density 
ratio in the nuclear region of NGC 4418 is a; 1. 

The PIMS lines of H 3 + are not detected in NGC 4418, but 
the model for C core in Fig. [4] uses the same T ro i — 500 K as 
in Arp 220, and the predicted faint PIMS lines are nevertheless 
quite compatible with the observed spectrum. The model pre- 
dicts wing-like features at the wavelengths of the 9g - 9g and 
12j,— 12[ 2 ortho PIMS lines, which could be present in the spec- 
trum at very low (but not significant) signal-to-noise. Therefore, 
the non-detection of the PIMS lines does not imply that T mt is 
low in NGC 4418, but according to our model can also be due to 
moderate continuum opacity and dilution of the nuclear contin- 
uum emission at long far-IR wavelengths. 

3.4. Arp 220 
3.4.1. OH + 

In Arp 220, the high-lying OH + 4j <— 3/< lines shown in Fig. [2^ 
are tentatively associated with the luminous western nucleus, 
C wes t, but they can be generally interpreted as the joint emission 
from the two nuclei (Paper I). In our model for C wes t, we include 
a velocity gradient through the absorbing shell of 130 km s , to 
reproduce the observed line broadening. We use two approaches 
to bracket the parameters inferred in Paper I for C wes t- For a 
compact/warm region, with 7d us t = 110 K and d — 93 pc, the 
high-lying OH + 4j <— 3j< lines are well fitted with a column 
of Af(OH + ) = 1.1 X 10 16 cm 2 , similar to the value found in 
NGC 4418. The predicted spectrum is shown with blue curves 
in Fig. |2 The OH + absorption in the 4j <— 3y lines, however, 
does not require such high 7d us t- For a colder and more extended 
region, with 7d us t = 90 K and d = 160 pc, the column density 
is A^(OH + ) = 4.5 x 10 15 cm 2 . These models, however, do not 
reproduce the strong 2j <— 1 /< and 32 <— 22 lines (Fig. |2^-d), 
so another component is required in addition. We identify this 
additional component as C ex tended (Paper I), which has moder- 
ate excitation in H2O and OH. The C exte nded component can be 
identified with the (inner region of) the external disk traced by 
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Table 1. Line equivalent widths and absorption fluxes. 



Line 



NGC 4418 Arp 220 



OH + 2i - 1 
OH + 2 3 - 1 2 
OH + 2 2 - 1 2 
OH+ 2 2 - 1, 
OH+2, - 1, 
OH+ 4 5 - 3 4 
OH + 4 3 -3 2 
OH+ 4 4 - 3 3 
pH 2 + 2 12 |- loif 
P H,0 + 2 12 |- lot I 
P H 2 + 2 12 |- loif 
P H 2 + 2 21 |- 1 10 1 
pH 2 + 2 21 |- 1 10 | 
oH 2 O+3 I3 |- 2 02 f 
oH 2 O + 3 13 |- 2 02 | 
oH,0 + 322? -2u| 
oH 2 + 3 22 |-2„| 
0H7O 4o4 2 _ 3i3| 
oH 2 O + 4 04 |-3 13 | 
oH 3 + 3- - 3 + 
pH 3 + 4--4 4 + 
pH 3 + 5- - 5+ 
oH 3 + 6 6 - 
pH 3 + 7" - 7+ 
pH 3 + 8 8 - 8+ 
oH 3 + 9, - 9J 
P H 3 0+ 107 - 10+ 
oH 3 + 127, - 12+ 
P H 3 + 5- 4 - 4+ 
oH 3 + 4- - 3 + 
pH 3 + 4" - 3 + 
pH 3 + 4^-3| 
oH 3 O + 4 - 3 + 
pH 3 + 3- - 2\ 
pH 3 + 37 - 2\ 



F, 

Mower 




4 , 

/'e-l 'l 

l,s ) 


^eq 
^Km S ) 


w cm ) 


S ) 


( 1 0.-21 \\T ~ m 

\iv w cm 


43.6 


148.696 


0.10 


6.7(1.8) 


0.79(0.21) 


17.8(0.9) 


10.0(0.5) 


46.6 


152.989 


0.18 


b 


b 


b 


b 


46.6 


147.768 


0.048 


3.9(1.6) 


0.47(0.19) 


20.7(1.2) 


12.0(0.7) 


49.6 


152.369 


0.14 


19.5(2.1) 


2.1(0.2) 


45.3(1.3) 


24.6(0.7) 


49.6 


158.437 


0.071 


5.9(3.0) 


0.6(0.3) 


19.8(1.1) 


9.6(0.5) 


281.9 


76.510 


1.6 


3.1(0.7) 


1.6(0.4) 


5.8(0.4) 


10.6(0.7) 


282.5 


76.245 


1.5 


3.2(0.8) 


1.6(0.4) 


3.1(0.4) 


5.7(0.7) 


285.5 


76.399 


1 .5 


3.8(0.8) 


1.9(0.4) 


4.0(0.4) 


7.4(0.7) 


u.u 


1 O A 11(1 

184.339 


n not 


^ A O 

< 4.0 


< 0.25 


4/.0(Z.l) 


10.3(0.5) 


0. 1 


182.915 


0.073 


< 3.7 


< 0.21 


30.8(2.1) 


7.5(0.4) 


0.0 


182.640 


0.016 


< 3.7 


< 0.21 


14.5(1.6) 


3.6(0.4) 


29.1 


105.742 


0.46 


7.6(2.3) 


2.1(0.6) 


6.2(0.8) 


7.7(1.0) 


30.4 


104.719 


0.39 










89.1 


143.813 


0.18 


5.2(2.5) 


0.65(0.31) 


22.6(1.1) 


14.0(0.7) 


89.3 


143.270 


0.17 


< 3.0 


<0.38 


8.3(0.9) 


5.2(0.6) 


125.0 


89.590 


0.59 


7.2(2.3) 


2.9(0.9) 


10.2(0.7) 


16.6(1.1) 


125.9 


88.978 


0.57 


5.4(1.6) 


2.2(0.6) 


2.2(1.1) 


3.6(1.8) 


189.2 


145.917 


0.16 


4.4(0.8) 


0.55(0.10) 


2.0(0.5) 


1.2(0.3) 


189.8 


146.217 


0.15 


4.7(0.8) 


0.58(0.10) 


5.5(0.6) 


3.2(0.4) 


102.0 


180.209 


0.083 


< 3.9 


< 0.25 


28(3) c 


8.5(1.0) c 


179.0 


178.994 


0.091 


< 3.9 


< 0.27 


12.7(1.3) 


4.0(0.4) 


273.4 


177.272 


0.097 


< 6.0 


< 0.44 


13.9(1.6) 


4.4(0.5) 


385.1 


175.063 


0.10 


< 4.6 


<0.38 


18.6(1.3) 


6.4(0.4) 


514.0 


172.396 


0.11 


< 3.4 


<0.28 


11.8(1.4) 


4.2(0.5) 


660.1 


169.303 


0.12 


d 


d 


13.7(1.2) 


5.4(0.5) 


823.3 


165.819 


0.13 


8.0(2.0) 


0.71(0.18) 


16.5(2.0) e 


6.9(0.8) c 


1003.5 


161.984 


0.14 


< 3.7 


<0.36 


7.5(1.1) 


3.3(0.5) 


1414.7 


153.424 


0.17 






8.5(1.5) f 


4.5(0.8) f 


179.0 


60.164 


0.49 


5.8(2.0) 


4.1(1.4) 


5.5(1.2) 


12.1(2.7) 


102.0 


69.538 


0.38 


6.1(1.3) 


3.7(0.8) 


17.0(1.5) 


33.6(2.9) 


139.3 


70.254 


0.63 


< 2.3 


< 1.4 


5.4(0.7) 


10.4(1.4) 


161.7 


70.684 


0.77 


5.9(1.9) 


3.5(1.1) 


4.8(0.5) 


9.1(1.0) 


169.1 


70.827 


0.81 


8.0(2.0) 


4.6(1.2) 


7.7(0.8) 


15.8(1.5) 


42.3 


82.274 


0.28 


4.0(1.3) 


1.8(0.6) 


7.9(0.6) 


13.8(1.0) 


64.7 


82.868 


0.44 


< 1.2 


< 0.6 


13.1(0.7) 


23.2(1.2) 



Numbers in parenthesis indicate lcr uncertainties. 
The line is strongly contaminated. 
Contaminated by H 2 and CH (Fig.gji). 
Dominated by HCN (Paper I). 

Contaminated by NH 3 ; the line flux in Arp 220 is estimaded by subtracting the modeled NH 3 flux from the flux of the whole feature (Fig.|4g). 
Tentative detection, at the wing of the NH 153 jim line (Figs. [2^ and[4j). 



CO dDownes & Eckartll2007t) . In Paper I we estimated d ~ 650 
pc, though this is uncertain by at least a factor of 2. The addi- 
tional OH + absorption could indeed be produced in a relatively 
small region (d < 300 pc) around the nuclei; our results are 
relatively uncertain for this component. Assuming that the ex- 
tended OH + mainly absorbs the continuum of the C ex tended com- 
ponent (Fig. 1 in Paper I), it most probably has a column of 
N(OH + ) * (2 - 4) x 10 15 cm" 2 . This reflects the widespread 
presence of OH + in the regions surrounding the nuclei. 



The joint contribution of C west and C ex tended to the OH + ab- 
sorption, shown with a red curve in Fig. [2] yields a reasonable 
fit to most OH + lines in Arp 220, though the 32 <— 22 line at 
104 fim remains underpredicted (and blended with an H2O line). 
The wing-like spectral feature at 153.45 //m is not reproduced, 
though it is probably due to H 3 + (Sect. 13.4.31 ). 



3.4.2. H 2 + 

The above models for C wes t were applied to H20 + . Using the 
parameters for the compact/warm approach, Ti mt = 110 K and 
d = 93 pc, we derive a H20 + column density similar to that 
found for OH + , N{H 2 + ) ~ 0.9 x 10 16 cm 2 . The model result is 
shown in Fig.[3](blue curves). In this model the 322211 5/2 <— 
3/2 and possibly the 331220 7/2 <— 5/2 lines are overpredicted, 
but the other ortho-H20 + lines are well reproduced. Using the 
more extended, colder approach, with T^st = 90 K and d — 160 
pc, the column density is Af(H 2 + ) = 2.5 x 10 15 cm" 2 . In the 
nuclear region of Arp 220, the most likely OH + -to-H20 + ratio is 
1 - 2. 

Apart from the para-H20 + ground transition, there is only 
one excited, non-contaminated para-H20 + line detected in the 
spectrum of Arp 220, the 221110 5/2 - 3/2 at 105.7 fjm. As 
can be seen in Fig. |3j), the absorbing flux in this line is well 
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Table 2. Inferred column densities, OH + abundances, and OH + /H20 + and OH + fHj,0 + column density ratios. 



Source 


NGC 4418 


Arp 220 


Component 


C -C 

y - core warm 


Cextended 


c 

'-'west 


Cextended 


dust (K) 


110- 150 


variable 1 -' 


90- 110 


variable 


N(OH + f (10 16 cm" 2 ) 


0.6- 1 


~ 0.3 


0.5- 1.1 


2-4 


iV(H 2 + ) a (10 16 cm" 2 ) 


0.3 - 0.4 




0.25 - 0.9 


~0.5 d 


N(H^O + f (10 16 cm" 2 ) 


0.5 - 0.8 




0.9 - 2.7 e 




^(OH + ) b (10- 8 ) 


1.5-2.5 


~ 2 


1.1 - 2.8 


~ 4 


OH + /H 2 + 


1.5-2.5 




1.0-2.0 


5-10 


OH + /H,0 + 


~ 1 




~ l c 





a Columns toward the nuclei are calculated within a screen approach. 

b Estimated OH + abundance relative to H nuclei, using a normalization column density of A'(H) = 4 x 10 23 cm 2 for the nuclear components 
(Paper I). 

c ?dust varies with radial position as it is calculated from the balance between heating and cooling. 

d The p-H 2 + ground-state lines at 183 - 184 yum in Arp 220 require an additional absorbing halo component (C na io) with rV(H 2 + ) ~ (2.5 - 
3.0) x 10 15 cm- 2 . 

c In Arp 220, the total H30 + column density is twice the nuclear OH + column, but since the H30 + lines are redshifted relative to the OH + lines, 
we roughly estimate an OH + /H 3 + ratio of 1 in the nuclear region where the OH + lines are formed. In the C ex tended component, the OH + /H 3 + 
ratio is uncertain. 



reproduced with our models, which use an ortho-to-para ratio 
of 3, the high-temperature limit. This provides a slightly better 
fit than the one obt ained with an ortho- to-para ratio of « 4.8, 
the value derived bv lSchilke et al.l(l2010h toward Sgr B2, which 
underestimates the absorbing line flux in Arp 220 by 35%. 

The ground-state p-H 2 + 212101 lines at 183 - 184 yum can- 
not be reproduced with the C west component (Fig. [3^). Models 
for the Cextended were also unsuccessful in reproducing these 
spectral features, as the required high (1.6 x 10 16 cm 2 ) H20 + 
columns have also the effect of predicting too strong absorption 
in other lines, specifically in the 322211 components. The p- 
H20 + 183-184 fim lines peak at a velocity different from the 
H20 + HE lines (Fig. [6), and are thus likely generated in a low 
excitation region that generates faint absorption in other lines. 
This component is likely to be the same as that producing strong 
absorption in the ground-state lines of H2O and OH, i.e. the Chaio 
component (Paper I). The spatial extent of this component is not 
constrained, and can be interpreted as foreground gas located 
along the line of sight to the nuclei absorbing their nuclear far-IR 
emission. As shown below (Sect. [3~4~3l ). however, it probably in- 
cludes gas within a few hundreds pc from the nuclei; the C ex tended 
and Chaio components are hard to separate physicalljQ. The de- 
rived H20 + column in this component is A^(H20 + ) * 3 x 10 15 
cirr 2 , and the fit is shown with a light blue curve in Fig.[3h. 

The amount of H20 + in C exter ,ded is better traced with the ex- 
cited H20 + detected in emission with Herschel/SPIRE by Rll. 
In Fig. [8] we show the predicted SPIRE spectrum of H20 + in 
Arp 220 between 600 and 1200 GHz. The contributions by 
the C W est, Cextended, and Chaio are shown with blue, green, and 
light-blue curves (vertically shifted for clarity), respectively, and 
the red curve shows the total absorption/emission. At these fre- 
quencies, the contribution by the nuclear region (C wes t) is not 
expected to be dominant. The absorption in the ground-state 

8 The difference between the C ex tended and Chaio components in our 
models for Arp 220 arises from radiative transfer effects. The C ext ended 
is characterized by optically thin, extended far-IR emission, where 
the molecules are still significantly excited via absorption of locally- 
emitted far-IR photons. The Chaio, however, describes foreground, low- 
excitation gas with strong, optically thick emission behind arising from 
the nuclei. The main difference may simply be the location of the 
molecules relative to the far-IR optically thick nuclei and the observer, 
thus basically describing the same spatially extended gas. 
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Fig. 8. Predicted SPIRE spectrum of H 2 + in Arp 220 between 
600 and 1200 GHz. The contributions by the C wes t, C ex tended, and 
Chaio are shown with blue, green, and light-blue curves (verti- 
cally shifted for clarity), respectively, and the red curve shows 
the total absorption/emission. 



111000 / = 3/2 <- 1/2 and 1/2 ^ 1/2 o-H 2 + lines at 
1 100 - 1 150 GHz is expected to be dominated by the Chaio com- 
ponent, also responsible for the absorption in the ground-state p- 
H20 + lines at 183 - 184 /mi. Most relevant is the predicted emis- 
sion in some H20 + lines at 700 - 800 GHz, which are expected 
to arise from the C ex tended component. The strongest H20 + emis- 
sion lines are the 2o2 — > In 5/2 — > 3/2 at 742.1 GHz and the 
2n — > 2o25/2 — > 5/2 at 746.5 GHz, in agreement with the emis- 
sion features reported by Rl 10. Two other transitions with signif- 
icant predicted emission above the continuum are the 2o2 — > 1 1 1 
3/2 1/2 at 721.8 GHz and the 2 n 2 02 3/2 3/2 at 761.9 
GHz. It is worth noting that the modeled emission features are 
generated through radiative pumping and not through collisions 



9 Rl 1 assigned the 746.5 GHz feature to the 2 02 -> 1 1 1 3/2 -> 3/2 at 
nearly the same frequency, but this transition has an Einstein coefficient 
5.4 lower than the main 5/2 — » 3/2 component and is expected to give 
weak emission 
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-collisional excitation is not included in our calculations. With 
a column density in C ex tended of N(H20 + ) ~ 5 x 10 15 cm -2 , the 
lines fluxes at 700 - 800 GHz reported by Rll are nearly re- 
produced. Our models thus favor an OH + /H20 + column density 
ratio of > 5 for C exte nded, significantly higher than in the nuclear 
region (Table|2]). In the Chaio, the OH + /H20 + is uncertain as the 
OH + ground-state lines are certainly saturated (Rl 1). 

3.4.3. H 3 + 

The H30 + lines in Arp 220 show an intriguing behavior: on the 
one hand, non-metastable RI lines are detected (Fig.|4jmnp); on 
the other hand, although they are radiatively pumped by far- 
IR photons to moderate excitation, these lines peak at roughly 
the same velocity as the LE lines of H2O and H20 + transitions 
(Pig. O, i.e. they are redshifted relative to the HE lines of other 
species. 

From the population diagram in Fig. [5] we have used two 
values for T ml : 180 K up to the K = 5 ladder, and 500 K for 
higher-lying metastable levels. The blue lines in Fig. [4] show the 
model results with T iust = 90 K, d = 160 pc, and N(H 3 + ) = 
9 x 10 15 cm 2 . The model for H 3 + nearly accounts for the 
pure-inversion metastable lines (panels a-i), which have moder- 
ate opacities (r ~ 0.3-1), though the 3j <— Saline is slightly un- 
derpredicted. The rotation-inversion lines (panels j-p) are more 
optically thin. 

Figure [9^ shows the flux ratios of the pure inversion, 



metastable 6 fi 
to the 47 



6+ 8; 



9+ and 10 



10 



10| lines 



t j 1 4^ line (solid lines) as a function of T mt , and com- 
pares them with the measured values (dashed lines). The flux 
measured for the 9 9 <— 9g line has been estimated from the total 
flux of the broad 165.8 fim feature by subtracting the modeled 
contribution due to NH 3 (Fig. gg). The K = 8/K = 4 (blue) 
and K = 9/K = 4 (green) ratios indicate T mt ~ 600 - 800 K, 
while the K = 6/K = 4 (red) and K = IQ/K = 4 (black) fa- 
vor ~ 450 K. The high T rot inferred from the H 3 + lines is in 
striking contrast with that inferred from NH 3 (Fig. [9J5): the NH3 
line ratios indica te r rot ~ 150 K, in general agreement with the 
value inferred bv lOtt et aU d201 ll) from the pure inversion lines, 
and with the HCN excitation (Paper I). Either the two species 
are sampling different regions, and/or r, ot as derived from H 3 + 
better reflects the formation process rather than T gas . 

The strength of the non-metastable rotation-inversion lines is 
sensitive to 7d us t, and in Fig. [9]; we compare the fluxes of some 



of these lines relative to that of the metastable 3 



2 + 
-2 



line. We 



focus in this analysis on those lines that have similar line shapes, 
i.e. the 4q <— 3y,4~ <— 3|, and4j <— 3^ transitions. Comparison 
with the measured values favors r^st ~ 90 - 110 K. The most 
reliable comparison is between the 4^ <— 3J and the 3j <— 2J, 
as they belong to the same K-ladder, favoring T^ usl ~ 90 K. 
Furthermore, due to the warmer SED, Tjust = 110 K generally 
underestimates the PIMS lines relative to the RI lines. In con- 
trast, the H2O HE lines are compatible with higher Ta ust (Paper 
I), and detection of the NH3 6 <— 5t and 6^ <— 5| lines at 
83.85 fim (Paper I and confirmed with the new data set), with 
£iow ~ 400 K and pumped through radiation, better supports 
Tdust = HO- 130KforNH 3 . 

The value of Af(H 3 + ) depends on the size of the far-IR con- 
tinuum source. To produce the observed H 3 + absorptions, the 
models strongly favor optically thick far-IR continuum emission 
up to at least 180 /urn. An upper limit to the size would then be 
230 pc, as the continuum source with 7d ust = 90 K would pro- 
duce the entire bolometric luminosity of the galaxy, 1.5 X 10 12 
L . This upper limit is, however, unreliable, as C ex tended is ex- 
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Fig. 9. Model results for H 3 + and NH 3 in Arp 220. a) Ratios of 
pure inversion metastable lines of H 3 + versus T lot : K = 6/K = 
4 (red), K = 8/K = 4 (blue), K = 9/K = 4 (green), and K = 
4 (black), b) For NH 3 , the fluxes of the 6 5 - 5+ (red), 



10/isr 

7g - 6 6 (green), and 8 6 - 7g (blue) are normalized to the flux of 
the 124.9 fim feature (composed of the 4o - 3o and 4^ - 3| NH 3 
lines, Paper I), c) Ratios of rotation-inversion lines of H 3 + as a 
function of 7d us t- Dashed lines in all panels indicate the measured 
values, with typical uncertainties of 20%. 



pected to generate some far-IR emission, and there is also an 
important fraction of the luminosity associated with dust warmer 
than 90 K. Sizes smaller than 100 pc are also unlikely, as the 
H 3 + PIMS lines would become optically thick, and the fit 
to the whole H 3 + SLED would get worse (specifically the 
6 6 <— 6g line would be underpredicted due to saturation). We 
estimate d = 105 - 160 pc, overlapping with, but at the high 
end of our derived size range for the nuclear OH + and H20 + 
(93 - 160 pc) and significantly smaller than our estimate for the 
extended region as derived from the continuum emission (~ 650 
pc, Paper I). This size range implies a column density range of 
N(H 3 + ) = (2.7-0.9)xl0 16 cirT 2 . These columns areafactorof 
2 higher than those of OH + for the same source sizes. However, 
the H 3 + transitions do not have the same line shapes as the OH + 
lines (Fig.|6h). There is, nevertheless, a good match between the 
OH + HE and H 3 + 70.8 fim profiles on the blueshifted side, 
suggesting some spatial overlap between both components. On 
these grounds, we favor a column density ratio OH + /H 3 + ~ 1 
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in the nuclear region of Arp 220, with an uncertainty of ±50%. 
The relatively large size found for the H30 + absorption, the shift 
of the peak of the H30 + velocity profile towards those of the 
low excitation lines, and the moderate T^ usl found for the H30 + 
lines, suggest that the H30 + lines, while nuclear, trace a transi- 
tion region between the high-excitation (HE) component and the 
extended component. 

The model in Fig.|4]predicts a flux for the 3£ — > 2 2 365 GHz 
line of 1.1 x 10~ 23 W cirT 2 , a factor of 1.5 lower than measured 
by vdT08. The modeled line is slightly inverted but with negli- 
gible amplification, and thus behaves as optically thin with the 
flux proportional to the column. Either it is sampling in addition 
deep regions that are extincted in the far-IR, or it is sampling a 
larger region (vdT08), because the 365 GHz emission line does 
not require far-IR optically thick continuum emission behind. 

3.5. Summary 

The analysis of the O-bearing molecular ions in NGC 4418 
and Arp 220 reveal (f) high nuclear OH + column densities of 
(0.5 - 1) x 10 16 cirT 2 in both sources, with estimated abundances 
relative to H nuclei of (1 - 3) x 10~ 8 ; (if) OH + /H 2 + column 
density ratios of ~ 1-2.5, with NGC 4418 at the high end of the 
range; (Hi) high OH + columns of ~ (2.5 - 4) x 10 16 cirT 2 in the 
Cextended component of Arp 220, with an estimated abundance of 
~ 4 x 10~ 8 , as well as a high OH + /H 2 + ratio of ~ 5 - 10; (fv) 
H 3 + column densities of ~ (0.5-0.8)xl0 16 cirT 2 in NGC 4418 
and ~ (0.9 - 2.7) x 10 16 cm" 2 in Arp 220, the latter essentially 
nuclear but likely tracing a region slightly more extended than 
that traced by OH + . In the nuclear region where the OH + lines 
are formed, OH + /H30 + ~ 1 is estimated. 



4. Chemistry 

4.1. Framework 

The O-bearing molecular ions are expec ted to be formed in 
gas irradiated by cosmic rays or X-ra ys (Malonev et al., 1996; 
iGerin et al.L 120101; iNeufeld et all 1201 0t lMeiierinketaiIl2011l) 
The sequence is initiated by the ionization of H atoms and 
H 2 molecules; OH + and H20 + are subsequently produced by 
ion-neutral chemistry, and their abundances are sensitive to the 
gas ionization rate £ and the molecular fraction /h 2 . In re- 
gions with a significant atomic fraction, H + transfers its charge 
to atomic oxygen, and the + reacts with H 2 to form OH + ; 
H 2 + is the n formed from OH + + H 2 , and similarly H30 + from 
H 2 + + H 2 (IGerin et all I20I0L INeufeld et al., 2010). In regions 
with high molec ular fraction, OH + and H 2 Q + can also be formed 
via Hi + O re.g.)Herbst & KlemDeredfT973l: IGerin et all 120101: 

with 



IHollenbach et all 120121) . Furthermore, in environments 
high reservoirs of OH and H 2 0, OH + and H 2 + can be produced 
from the charge exchange reactions OH + H + an d H 2 Q + H + , 
as we ll as from photoionization of OH and H 2 dGupta et all 
2010), while H 3 + may also be formed from H 2 + H+. Finally, 
H + can also be produced through C + + OH and CO + + H, and 
this F UV/chemical route does not require X/cosmic ray irr adia- 
tion dSternberg & Dalgarnol 119951: IHollenbach et alll2012T) . 

In order to discriminate between the above paths of OH + , 
H 2 + , and H30 + production, and to infer the physical condi- 
tions implied by the observations, we h ave run a grid o f chem - 
ical models with the code described b y Bruder er et al. | (|2 009). 
based on previous works by Staub er et al.l (120051) and lDotv et al.1 
(2002). The chemical models calculate the steady state abun- 
dances of all relevant species based on the UMIST06 reaction 
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Fig. 10. Predicted (a) OH + /H 2 + and (b) OH + /H 3 + ratios as 
a function of Xe X 10 3 //h 2 - Symbols show results for all gen- 
erated che mical mode l s, whi le th e solid lines show th e analytic 
curves bv lGerin et all (1201 Oh and INeufeld et alj d2010h . The in- 
ferred OH + /H 2 + ratios in the nuclear regions of NGC 4418 and 
Arp 220, and in the C ex tended of Arp 220, are indicated with the 
shaded dark and light regions, respectively. Gas temperatures of 
150 K (black) and 550 K (red) are assumed. 



rates (IWoodall et all fcOO?). In our simple models, the gas with H 
nuclei density «h and molecular fraction /h, is directly exposed 
to a cosmic ray and/or X-ray flux that produces a total ioniza- 
tion rate per H nucleus £ (including secondary ionizations). We 
ignore any other external agent (e.g. dissociating UV radiation) 
in the calculations, though internally generated UV radiation is 
included. 

In their models for the prod uction of the O-b e aring molec- 
ular ions in interstellar clouds, IHollenbach et al ] (120 121) found 
that PAHs have important effects on their expected abundances. 
Specifically, a fraction of H + ions produced by cosmic rays was 
found to be neutralized mostly by PAH , truncating the path for 
the production of OH + , and thus reducing the effi ciency 6qh+ of 
OH + formation relative to previous calculations (INeufeld et all 
2010). Recent obs ervations of both O H + and HJ toward the 
Galactic W5 1 by llndriolo et al.l (1201 2|) have indeed indicated a 
relatively low value of £qh+. In NGC 44 18, however, PAH emis- 
sion is very weak dSpoon et all 120071) . and it is al so weak in 
Arp 2 20 relative to the submillimeter continuum dHaas et all 
120011) . Although this may be due to extinction of nuclear far-UV 
photons or foreground mid-IR extinction, rather than low PAH 
abundance, we conservatively switch-off the charge exchange 
between gas and dust grains in our calculations. In this respect, 
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the ionization rate £ we derive below (Sect l4.5b should be viewed 
as a lower limit. 

Under our sim plified assumpt i ons, and as shown by pre- 
vious studies (e.g. iMalonevetall 119961: iNeufeld et al.L HoTot 
iHollenbach et alll2012l) . the relevant input parameters are <57«h, 
fn 2 , the gas temperature (r gas ), and the oxygen abundance in the 
gas phase (xo)- For the latter parameters we adopt r gas = 150 
and 550 K, and solar metallicities with^o = 3 x 10~ 4 . Freeze- 
out of atomic oxygen on grains surfaces is thus ignored owing 
to the high Tanst in the nuclear regions of these sources (e.g. 
IHollenbach et all 120121) . 

In the calculations, we treat /h, as an independent param- 
eter, even though its value can in principle be obtained from 
the H2 formation and destruction rates. Treating /h 2 as a free 
parameter has the advantage that it allows us to run simple 
single-point models without additional assumptions on the ge- 
ometry or H2 formation efficiencies. Furthermore, fn 2 also de- 
pends on other physical processes that are difficult to evaluate 
in the sources under study: photodissociation by UV radiation, 
and the evolutionary state of the cloud in low density (diffuse) 
cloud s, where /h, only attains steady state after > 10 Myr (iLiszti 
120071) . Concerning the latter point, the ground-state lines of OH + 
and H20 + have been detected in Galact i c diffuse clouds indi - 
cating a very low / H , (iGerin et all 1201 Ot INeufeld et all 1201 Ol) , 
probably reflecting regions that are evolving in /h 2 . Our calcu- 
lations are intended to describe the observations of both the nu- 
clear regions and the extended component in Arp 220; in the 
latter case similar conditions to those found in Galactic diffuse 
clouds may be present. For a given £/«h, the maximum value 
that /h, can attain is described ap 



Jh 2 



('♦£)' 



1 + 



3.7 x 10- 17 cm 3 s- 



(2) 



where we have adopted for the rate coefficient of H2 forma 
tion on grain surfaces the reference value 



Rc 



3 x 10~ 18 cm 3 s~ l and T. 



3 o-l 



gas 



R g t1 12 with 
150 K, and applied this 
constant value to all cases ( i.e. regardless of T m * and 7d U st) for 
simplicity (see discussion in Kaufman et al., 1999). For the £/«h 
values considered in the following sections, 5 x 10~ 19 , 5 x 10 -18 , 
2 x 10~ 17 , 5 x 10~ 17 , and 5 x 10~ 16 cm 3 s -1 , /™ x is estimated as 
0.99, 0.88, 0.65, 0.43, and 0.07, respectively. Values of f Hl be- 
low /™ x are obtained in our models by manually increasing the 
H2 destruction rate. Equation [2] appears to be roughly consistent 
with the models by IBavet et alJ yOll), who obtain ~ 0.7 
at high Ay for £/«h = 10~ 17 cm 3 s _1 and solar metallicities. 
However, the actual /™ x values could be lower, owing to the 
de creasing H? recombination efficiency in warm grains (> 100 
K; ICazaux & Tielensl 12004 ICuppen et all 120101) . The grid of 
models cover the f//iH - /h, plane, with /h 2 varying from 0.02 
to /™\ and £/n H up to 5 x 1 0~ 16 cm 3 s" 1 . 



4.2. The OH + /H 2 + and OH + /H 3 + ratios 

Figure[10]shows with squares, for all generated models, the pre- 
dicted Ri = OH + /H 2 + and R 2 = OH + /H 3 + ratios as a func- 
tion ofxe x 10 3 / /h 2 , wher e y e is the electron abun dance relative 
to H n uclei. As shown bv lGerin et alJ (1201 Ol) and INeufeld et alJ 
d2o!oh . if the H20 + is exclusively formed by reaction of OH + 



10 It is assumed that each H2 ionization gives rise to one net H2 de 
struction due to the expected dominant paths HJ + H 2 



H+ +e 



Ho + H. 



H+ + H and 



with H2, ^1 only depends onxd '/h 2 and r gas , regardless of the 
formation mechanism of OH + , according to: 



R x = OH + /H 2 Q + = 0.63 + 0.85 x 



X 10 3 



(3) 



where 7300 = r gas /300K. If in addition H30 + is exclusively 
formed from H20 + + H2, R2 is given by 



0.86x ^ 03 x(l + 1.34x ^° 3 V(4) 



R 2 = OH + /H 3 + 



These analytical dependences are shown with solid curves for 
T gas = 150 and 550 K in Fig. [TO] Eq. ® shows that a lower 
limit of R = 0.63 is expected when recombination of H20 + 
is ignorable and the unique route of H20 + destruction is re- 
action with H2 yielding H 3 + . For Xelfch ^ 10 -3 , recombi- 
nation of H20 + becomes increasingly important and R\ rises 
linearly with Xelfvh- This analytic approach to R\ is appropi- 
ate for regions where the formation of H20 + is dominated by 
H2 + OH + . However, at low^ e //H, values, i.e. in shielded/high 
density regions with low electron fraction and high molecular 
fraction, R\ deviates from the analytic curve because there are 
competing paths for H20 + formation (mostly H + + H2O and 
Ho + O) that do not involve formation of OH + , thus decreasing 
the OH + /H 2 + ratio below 0.63 (squares in Fig. [TOTt), On the 
other hand, R 2 varies strongly with Xelfu^ and significant de- 
partures from the analytical approach in eq. (0]l are also obtained 
for low^J/e,, and high r gas . 

Even in the nuclear regions of NGC 44 1 8 and Arp 220 where 
densities are high and high concentration of molecular species 
are observed (Paper I), we do find relatively high values of 
R\ and R 2 , though significantly lower than in Galactic diffuse 
clouds. R\ is found ~ 1 - 2 and could be even higher in NGC 
4418, and R 2 ~ 1 (Table|2]l. Thus the inferred ratios appear to 
indicate that the primary route for H20 + and H 3 + formation 
in the nuclear regions of NGC 4418 and Arp 220 is the same 
as in Galactic translucent clouds, H2 + OH + and H2 + H20 + , 
respectively, additional routes for H20 + and H 3 + formation 
being less important. We still expect, however, /h, to be sig- 
nificantly higher than in translucent clouds. Adopting conser- 
vatively /h, > 0.5, we obtain j e > 2 x 10~ 4 . This value is 
higher than that expected from complete ionization of carbon, 
1.4 x 10~ 4 , pointing to a significant, additional source of free 
electrons arising from H and H2 ionization. In the more extended 
and diffuse C ex tended of Arp 220, the high value of R[ probably 
reflects a relatively low value of /h 2 , resembling the situation 
found in Galactic translucent clouds. 

The full chemical results for OH + , H 2 0\ and H 3 + , to- 
gether with the corresponding abundances of electrons and H + , 
are shown in Appendix ICl where some analytical approaches to 
the abundance of OH + are also derived. We find that the most 
likely path for OH + formation is + + H2, with the contribution 
from the charge exchange reaction H + + OH and from direct 
OH ionization being of secondary importance even in the nu- 
clear regions, where high OH abundances are inferred (Paper I). 
Therefore, the sequence H + -> + -> OH + -> H 2 + -> H 3 + 
is favored as the production of the O-bearing cations. In 
Appendix [C] we also favor deep penetration of X-cosmic rays 
as the main source of H + production, with the contribution 
due t o the C + + OH -> CO| + H -> CO + H + FUV /chemical 
route dSternberg & Dalgarnol 1 19951: IHollenbach et all (20121) be- 
ing presumably important only in small volumes around O stars. 
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Fig. 11. Predicted OH + , H 2 + , and H 3 + abundances (left panels) and OH + /H 2 + and OH + /H 3 + abundance ratios (right panels) 
as a function of f /«h- In all these models, the molecular fraction is given by the maximum value /™ x in eq. ©. Gas temperatures 
of 150 K (black) and 550 K (red) are assumed. The ranges derived from the observations and models for the nuclear regions of 
NGC 4418 and Arp 220 are indicated by the shaded regions. 



4.3. The value of£/nn in the extended region of Arp 220 

In the Cextended component of Arp 220, where high OH + column 
and Ri values are found, our best fit values for £/«h and /h, are 
(0.5-2)xl0" 17 cm 3 s" 1 and 0.05 -0.25, respectively (Fig.lCTTV 
With this low molecular fraction, however, xoh+ is not far from 
saturation (see Appendix 0, which occurs at xoh+ Z 10~ 7 
and thus may be reached in some regions of C ex tended- The low 
value of /h, indicates that either the H 2 abundance has not fully 
evolved to steady state, or H 2 photodissociation keeps a large 
hydrogen fraction in atomic form. The electron abundance is ex- 
pected to be in the range Xe ~ (3 — 10) X 10~ 4 . 

These physical conditions are derived from the apparently 
high OH + column in C exte nded, and from the relatively low H 2 + 
column with the following caveat. The columns are derived from 
far-IR OH + absorption and submillimeter H 2 + emission lines, 
which may probe different regions. The lowest-lying absorption 
lines of H 2 + at 183-184 //m, observed in absorption toward the 
nuclei (ChaioX indicate a foreground H 2 + column of ~ 3 x 10 
cm" 2 . A fraction of the H 3 + observed toward the nuclei is ex- 
pected to arise also from foreground gas, with a column of also 
several x 10 15 cm 2 , though detailed models for PDRs irradi- 
ated by cosmic rays show that H 3 Q + attains significant col umns 
deeper into the cloud than OH + dHollenbach et all 120121) . The 



extended gas in Arp 220 likely harbors regions with a diversity 
of physical conditions. 

4.4. The value of£/n H in the nuclear regions 

In the nuclear regions of NGC 4418 and Arp 220, where ^oh+ ~ 
2 x 10~ 8 and R = 1-2.5 (Table 13, xor* is unsaturated 
(see Appendix ICt. i.e. results are sensitive to £/«h- However, 
the inferred £/«h depends on the assumed /h,. In these re- 
gions, we may expect steady state for H 2 , and we adopt /h 2 as 
given by eq. [2] (i.e. we assume that cosmic/X rays are the only 
source of gas dissociation/ionization). Figure Q~T] shows the cor- 
responding modeled cation abundances an d ratios. Model results 
can be summariz e d as follows (see also |N eufeld et al., 2010; 
iGerin et all l2010t iGupta et all l2010t iHollenbach et all 120121) . 
For very high £/«h, OH + can still attain high abundances even 
with low /h, < 10%, but H 2 + is primarily destroyed through 
recombination and the abundances of H 2 + and H 3 + are low. 
For low £/«h, low abundances of OH + and H 2 + are obtained, 
but the H 3 + abundance for r gas = 550 K remains high as H 3 + 
is formed through reaction of highly abundant H 2 with and 
H + . Abundances of OH + and H 2 + around 10~ 8 are obtained at 
typical £/«h ~ 10~ 17 cm 3 s _1 . The abundance of H 3 + is very 
sensitive to r gas , and 7? 2 is very sensitive to both r gas and f /«h- 
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The inferred abundances and ratios, marked by shaded re- 
gions in Fig. [TTJ are consistent with f/«n in the range (1 — 
2) x 1CT 17 cm 3 s _1 . The electron abundance is expected to be 
in the range Xe ~ (1-5 - 5) x 10~ 4 . These results are sensitive 
to the assumed H2 fraction, with higher (lower) /h 2 implying 
higher (lower) £/«h- Decreasing yn, in eq- © to 10~ 17 cm 3 
s _1 , the observed abundances and ratios can be explained with 
£/«h ~ 5 x 10~ 18 cm 3 s _1 . Further study of the H2 reformation 
on warm dust grains in these regions is required to refine our 
estimates. 

4.5. The total ionization rate per H nucleus (() 

Estimation of the total ionization rate per H nucleus ( requires 
knowledge of the gas density. In the C ex tended of Arp 220, an 
average density of «h ~ 10 3 cm -3 is expected from the contin- 
uum emission and estimated size (Paper I). Therefore £ ~ 10~ 14 
s _1 , about 100 times higher than in Galactic diffuse clouds. This 
value is however uncertain, as the absorption in the OFF lines 
may be selectively produced in regions with densities much 
lower than the estimated average, and thus f significantly lower 
than our above estimate is not ruled out. 

In the nuclear regions of NGC 4418 and Arp 220, very high 
densities are expected. In the C core of NGC 4418, densities of 
3 x 10 6 cm -3 are derived from the faint OH absorption in the 
ground-state lines at central velocities, from the pattern of H2O 
absorption, and from the high excitation of HCN (Paper I). If 
the observed OH + were formed in this region, this would imply 
an extreme, very unlikely value of f ~ 10~ 10 s . However, it is 
plausible that the OH + lines are formed in a more tenous, nuclear 
region. There are several observations indicating that the molec- 
ular ions are formed in a component different from that respon- 
sible fo r the emission/ab sorption observed in neutral molecular 
species. iGupta et aH (l2010h have detected broad blueshifted ab- 
sorption in the ground-state lines of OH + and H20 + with HIFI 
toward the Orion KL, and their preliminary models appear to 
point toward a component of the low-velocity Orion outflow 
that has a density (~ 10 3 cm 3 ) far below the den s ities in ferred 
from other species. In diffuse gas, iNeufeld et al.l (l2010h show 
that the velocity distribution of the OH + and H20 + absorption 
is dissimilar from that of H2O; the absorption in the molecular 
ions (and atomic H) being more broadly distribute d in v eloc- 
ity space. Theoreticall y, the models by Ba vet et al] d201 lh and 
Meii erink et al.l d201 lh indicate that the OH + abundance is anti- 
correlated with the abundance of NH3 and HCN. Furthermore, 
the models by iMeiierink et alj d201 lb indicate that the H 2 
abundance is antic orrelated with £7«h> an d thus with the abun- 
dances of OH + and H20 + as well. 

We thus conservatively propose that the absorption observed 
in the high-lying lines of OH + and H20 + are formed in a rel- 
atively low-density region different from the dense component 
that accounts for the absorption in the high-lying lines of H2O, 
HCN, and NH3 (Paper I). This low density component could 
be interpreted in terms of an nuclear interclump medium, or a 
shell surrounding the high density nuclear gas. However, this 
component cannot be very extended as otherwise the far-IR ra- 
diation density could not excite the molecular ions up to the 
observed rotational levels (Slower ~ 285 K). The width of the 
OH + shell should be small in comparison with the radius of the 
far-IR source, which provides a lower limit on the density. Our 
models favor densities of at least ~ 10 4 crrT 3 in Arp 220, and 
~ 2 x 10 4 cm~ 3 in NGC 4418. For a density threshold of 10 4 
cm 4 , £ > (1 -2) x 10" 13 s 4 . In NGC 4418, where high densities 
were inferred in Paper I, this is a strong lower limit, and yet still 



f > several x 10 2 times the hi ghest values estimated in the Milky 
Way (severalxlQ- 16 s" 1 ; e.g.llndriolo eTdll2007tlNeufeld et all 
l2010HHollenbach et alll2012l:llndriolo et all 1201 2b . " 

5. Discussion 

5.1. The O-bearing cations and the neutral species 

Under the simple assumptions used in our schematic chemical 
models, predictions for the abundances of the neutral species 
OH, H2O, NH3, and HCN analyzed in Paper I are also obtained. 
The main results for the neutral species are: (z) OH attains abun- 
dances > 10~ 6 even for the lowest ^/«h = 5x 10~ 19 cm 3 s _1 and 
a wide range in fi{ 2 , in rough agreement with the derived values. 
(ii) H2O only attains the inferred high abundances of ~ 10~ 5 
for high Tgas and ^/«h < 2 x 10~ 17 cm 3 s _1 ; otherwise gas- 
grain chemistry (not included in our simple chemical calcula- 
tions) is required to account for the high H2O abundance (Paper 
I). (Hi) HCN only approaches abundances as high as ~ 10~ 6 for 
the lowest £/«h = 5 x 10 -19 cm 3 s _1 and high r gas , suggest- 
ing that indeed it is formed in regions with densities higher, and 
probably more protected from the ionization source, than those 
responsible for the molecular ions, (z'v) NH3 attains abundances 
of only 10~ 8 , far below the inferred ones, further indicating the 
importance of gas-grain chemistry. In general, clumpy regions 
with high contrast in nn, fii 2 , and £, and gas-grain chemistry, are 
most likely required to account for the full molecular absorption 
detected in the nuclear regions of both galaxies. The neutral/ion 
connection will be further analyzed in a future study. 

5.2. X-rays or cosmic rays? 

In the starburst galaxies NGC 253 and M82, detection of y- rays 
dAcero et all 120091; lAcciari et all 120091; lAbdo et all 1201 Oh in- 
dicate a high CR energy density of ~ 200 - 300 times the 
value (~ 1 eV/cm 3 ) in the Milky Way ( Abramow ski et alll2012l: 
iPersic & Rephaelill2012h . Since the y-rays are thought to be the 
result of cosmic rays accelerated by supernova remnant shocks 
that interact with high density target material, a high rate of 
cosmic -ray ionization indicates a high rate of supernovae in 
the central region of starburst gal axies (vsn ~ 0.03 yr _1 in 
NGC 253, lEngelbracht et all 1 19981) . which would be expected 
to be accompanied by observable galactic o utflows. These are 
indeed seen in both NGC 253 and M82 (e.g ISugai et all 120031: 
IStrickland & HeckmanLl2007tlVeilleux etalll2009h . " 

The cosmic -ray energy dens ity inferred in Arp 220 from both 
the supernova rate (~ 4 yr~\ Lonsd ale et all [2006) and syn- 
chrotr on radio emission i s ~ 5 00 - 1000 times the Galactic 
value (IPersic & RephaeliL 1201 2h . in good agreement with our 
lower limit of several x 10 2 times the Galactic value for The 
blueshifted component in Arp 220 traced by the high-excitation 
lines of H2O, OH, NH3, and HCN may indicate outflowing gas 
with a line of sight velocity of ~ 75 kms~' (Sect. l2~4b . With an 
extent of ~ 80 pc, this would imply a (local) very high mass out- 
flow rate of ~ 900 x («h/10 5 cm 4 ) M yr _1 . However, the out- 
flowing gas is difficult to distinguish from the rotating gas for the 
relatively low outflow velocities in Arp 220 (Sect. 12.41 ). in stark 
contrast with the ~ 1 000 fans" 1 outflow velocities traced b y 
OH in other ULIRGs dFischer et all 120101: ISturm et all 1201 lh . 
Nevertheless, clear indication of outflowing gas is seen in the 
P-Cygni HCO + profile s observed at submillimeter wavelengths 
dSakamoto et all l2009h . and is also suggested at redshifted ve- 
locities in the [O 1] 145 //m line and in the low-lying OH and 
H2O emission component profiles (Paper I). 
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The above independent estimates of the CR energy densi- 
ties in Arp 220, and the signatures of outflowing molecular gas 
with low-moderate velocities, appear to be compatible with the 
scenario of a high de nsity of cosmic rays producing the ob- 
served molecular ions. Papadopoulos (2010) has proposed that 
very high cosmic ray energy densities, of the order of those we 
derive here, can turn the bulk of the nuclear molecular gas in 
ULIRGs into giant CR-dominated Regions (CRDRs) rather than 
ensembles of Photon-dominated Regions (PDRs). In this frame- 
work, the O-bearing cations can be considered excellent tracers 
of these regions. 

The case of NGC 4 418 is unclear. This highly buried 
dSpoon et al.L 1200 ll 120071) and highly excited (Paper I) source, 
where the bulk of the luminosity arises from a very small vol- 
ume, is deficient in synchrotron emissio n with respect to the 
far-infrared emission dRoussel et alll2003l) . Furthermore, we did 
not find spectroscopic signatures of outflowing gas; on the con- 
trary the far-IR lines observed with PACS are narrow (Paper I). 
Therefore, cosmic rays may only be of secondary importance. 
An expla nation of the s ynchr otron deficit based on a nascent 
starburst (Rou ssel et all 120031) can hardly account for the ex- 
cited OH + detected in the source, though observations of excited 
OH + in Galactic and other extragalactic sources are needed to 
check this point. Preliminarily, NGC 4418 is better interpreted 
as a buried AGN, and the molecular ions as produced in a X- 
ray Dominated Region (XDR). But given that Arp 220 can be 
considered as a scaled-up version of NGC 4418, and that the 
columns of the cations are found similar in the nuclear regions 
of both sources, the possible XDR nature of Arp 220 is not ruled 
out. 

While both NGC 4418 and Arp 220 are underlumi- 
nous in hard (2 - 10 keV) X-rays for th eir computed SFRs 
dMaiolino et all 120031; llwasawa et al.U2005l) . the X-ray absolute 
fluxes may suffice to produce the inferred ionization rate. The 
distance-corrected X-ray fluxes of 0.1 and 0.5 erg s _1 crrT 2 esti- 
mated at the surface of the nuclear regions in NGC 4418 and 
Arp 220 (Paper I), respectively, are lower limits due to fore- 
ground absorption in both galaxies (with A^h ~ 10 23 cm" 2 , 
Paper I) and by the molecular shell itself where the ions reside. 
Depending on the X-ray spectrum, local (i.e. unattenuated) X- 
ray fluxes of > 0.5 erg s _1 crrT 2 may well produce £ in excess of 
10 ! s '. 

The strong hig h-ionization Fe K line in Arp 220 
( llwasawa et all 120051) appears to indicate that the faint X-ray 
emission is not primarily due to high-mass X-ray binaries (as 
in starburst galaxies), but could arise either from an inter- 
nally shocked hot bubble of starburst origin (the pressure of 
which would produce the observed molecular o utflow), or from 
a dee ply embedded supermassive black hole (Iwa sawa et all 
2009). In the latter case, the molecular ions would then be 
formed in (somewhat extended) regions that are illuminated by 
X-rays from the AGN that can escape to the outer (nuclear) re- 
gions primarily through inhomogeneities in the nuclear media. 
The drawback of this alternative is that the far-IR molecular lines 
are observed in absorption against dust emission that is optically 
thick even at far-IR wavelengths, and thus the X-ray emission 
from a putative buried AGN would then be strongly attenuated 
in arriving at such nuclear, but "surface" regions unless some 
special geometry is invokecQ- Further theoretical and observa- 
tional studies are required to discriminate between CRDRs and 
XDRs. 



11 The extinction in the X-ray regime for E x < 10 keV is higher tha n 
at 200 /um for a gas-to-dust ratio of 100 by mass (e.g. lMontmerlell200lh . 



5.3. H 3 + : formation pumping or hot gas? 

We show in Appendix IB1 that the H30 + population diagram for 
the PIMS lines in Arp 220 can be understood in terms of forma- 
tion pumping described by rf orm ~ 1000 K, but that plausible 
collisional relaxation and the lack of knowledge of collisional 
rates make the c ase uncertain. Interesti ngly, detailed chemical 
models of PDRs dHollenbach et all 1201 2l) predict that H 3 + is 
primarily formed in fully molecular regions, while both OH + and 
H20 + peak closer to the PDR surface where /h 2 is significantly 
lower than unity. On the contrary, in Arp 220, the H?0 + PIMS 
lines peak at the redshifted velocity (LE) component (Fig. [6j, 
which is expected to be more widespread and less dense than the 

blueshifted velocity (HE) component. 

Nevertheless, it is worth noting that Rangwal a et all 
have detected high-/ rotational lines of CO, which the authors 
interpreted as arising from hot (~ 10 3 K), low density (~ 10 3 
cm" 3 ) gas that is mechanically heated. If the H30 + lines are 
formed in the same region as the CO lines, the presence of such 
a hot molecular gas layer in Arp 220 may indeed suggest that 
both formation pumping and high r gas are responsible for the 
observed PIMS absorption. Due to the low gas density asso- 
ciated with the redshifted (LE) component, inefficient cooling 
may keep the gas at high r gas , which has the additional effect 
of boosting the H30 + abundance (Fig ITTb . On the other hand, 
the nuclear gas associated with the blueshifted ve locity (HE^ 
comp onent is dense, and gas-grain thermal coupling (Goldsmith, 
2001) may be efficient in cooling the gas to the inferred 150-200 
K (based on NH3 lines, Fig. [9J)). Collisional relaxation of the 
H30 + metastable populations may proceed efficiently at this r gas 

(Fig. ED. 

For a density of 10 4 cirT 3 and ( = 10~ 13 s _1 , the ex- 
pe cted r gas is ~ 300 K and rises very steeply with increasing 
£ dBavet et all 1201 lb : this result suggests X/cosmic ray heating 
may be important, together with mechanical heating, in explain- 
ing r gas in the nuclear regions of NGC 4418 and Arp 220. 

6. Conclusions 

The observational results on OH + , H20 + , and H30 + in 
NGC 4418 and Arp 220 are: 

- Absorption lines from excited levels of OH + with Slower up 
to 285 K are detected in NGC 4418 and Arp 220. While the 
highest-lying lines have similar strengths in both sources, the 
low-lying lines are more prominent in Arp 220. 

- Lines of H20 + are also well detected in absorption in 
Arp 220 up to a lower level energy of ~ 200 K, with strong 
absorption in the ground-state p-H20 + transition at 183-184 
jjm. In contrast, in NGC 4418 the ground-state line is not de- 
tected and only several of the excited lines are detected with 
relatively low signal-to-noise. 

- The pure-inversion, metastable lines of H30 + are detected in 
Arp 220 up to a lower level of ~ 1400 K, indicating a rota- 
tional temperature of ~ 550 K. These lines are not detected 
in NGC 4418. However, the rotation-inversion lines, some 
of them arising from non-metastable levels, are detected in 
both sources. 

- We identify two main velocity components in Arp 220 sep- 
arated by ss 75 kms -1 . The high-lying lines of all neutral 
(H 2 0, OH, NH 3 , HCN) and ionized (OH\ H z O + ) molec- 
ular species peak at the blueshifted component, while the 
lowest-lying lines peak at the redshifted one. The excep- 
tion corresponds to the pure-inversion metastable lines of 
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H30 + , which also peak at the velocity of the redshifted 
component. The blueshifted line components may be trac- 
ing the outflow shown in the millimeter lines of HCO + by 
Saka moto et al.l d2009l) . and grain mantle evaporation, sput- 
tering, and neutral-neutral reactions at high r gas may all be 
contributing to enhance the gas-phase abundances of H2O 
and NH 3 . 

The observations are analyzed with radiative transfer model- 
ing. The conclusions stemming from the analysis are: 

- The excited lines of OH + and H20 + , and the rotation- 
inversion lines of H30 + are pumped through absorption 
of far-IR radiation emitted by dust, but the pure-inversion 
metastable lines of H30 + are either chemically and/or colli- 
sionally pumped. 

- OH + column densities of (0.5 - 1.1) x 10 16 crrT 2 are derived 
toward the nuclear region of both sources for spatial extents 
of 20 - 30 pc in NGC 4418, and 90- 160 pc in Arp 220. 
The estimated abundance is (1 - 3) x 10~ 8 in both sources. 
In addition, in Arp 220, absorption from an extended region 
is required to match the low-lying OH + lines, with an es- 
timated column of a few x 10 16 crrT 2 and an abundance of 
several x 10~ 8 . 

- The nuclear columns of H 2 + are (0.3 - 0.9) x 10 16 cnr 2 
in Arp 220, and (0.3 - 0.4) x 10 16 cm" 2 in NGC 4418, indi- 
cating OH + /H20 + ratios of 1 - 2.5 (the lower value applies 
to Arp 220). In Arp 220, the emission lines detected with 
SPIRE indicate an additional extended region with columns 
of several x 10 15 , similar to the column derived from the 
ground-state p-H20 + lines, indicating a OH + /H20 + ratio of 
5 - 10 in the extended region. 

- The H 3 + lines indicate columns of (0.5 - 0.8) x 10 16 cnr 2 
in NGC 4418, and (1 -2)x 10 16 cnr 2 in Arp 220. In Arp 220, 
the observed high excitation of the pure-inversion metastable 
lines may be due to both formation pumping and high r gas , 
and the H30 + lines are likely tracing a nuclear transition re- 
gion into the more extended, low excitation gas. It is sug- 
gested that collisional relaxation of the H30 + metastable 
populations is important in the dense nuclear gas of Arp 220. 
The H 3 + lines are redshifted relative to the high-lying lines 
of OH + with uncertain velocity overlap, and we estimate a 
OH + /H30 + ratio of ~ 1 in the nuclear region of both sources. 

Simple chemical models are used to interpret the results of 
radiative transfer modeling. Our main conclusions are: 

- The sequence H + + -> OH + H 2 + H 3 + is fa- 
vored as the primary route to explain the high columns of 
the O-bearing molecular ions toward the nuclear region of 
both sources. 

- The production of H + is most probably dominated by 
X/cosmic ray ionization, rather than by FUV/chemical 
means (Appendix [Cb. An ionization rate per hydrogen nu- 
clei density of £/«h ~ (1 — 2) X 10~ 17 cm 3 s _I is estimated in 
the nuclei of both galaxies. 

- The full set of observations and models lead us to propose 
that the molecular ions arise in a relatively low-density in- 
terclump medium but with density of at least 10 4 cirT 3 , in 
which case the ionization rate is £ > 10~ 13 s _I , a lower limit 
that is over two orders of magnitude greater than the highest 
rate estimates of several x 10~ 16 s _I for Galactic regions. 

- The cosmic-ray energy density inferred in Arp 220 from both 
the supernova rate and s ynchrotron radio emission is ~ 10 3 
times the Galactic value (Persic & Rep haeliL l2012h . which is 



compatible with the lower limit found for f from the present 
observations and analysis. Nevertheless, an important contri- 
bution due to X-ray ionization in Arp 220 is also plausible. In 
NGC 4418, an X-ray Dominated Region is favored. Further 
theoretical and observational studies are required to discrim- 
inate between X-ray and cosmic ray ionization. 
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Appendix A: Excitation of OH + and H 2 + through 
formation pumping 

Here we show that formation pumping is not important for the 
OH + and H20 + level population distribution or absorption line 
fluxes in Arp 220 and NGC 4418. The effect of formation pump- 
ing on the excitation of O H + can be estimated fro m the statistical 
equilibrium equations bv lBruderer et alJ (1201 Oal) (eqs. B1-B3 in 
their Appendix B), by ignoring the collisional and induced emis- 
sion terms: 



dn 



NJ 



dt 



- 2_j n N+\J'A-N+\J';N,J - npjj y | Ajyr JJST-IJ' 
J' J' 

+ r«y - n N jK N j = 0, 



(A.l) 



where Asj-^pj' is the Einstein coefficient for spontaneous emis- 
sion from level (N, J) to level (N', J'), and T^j (crrT 3 s ) and 
Anj (s ) are the chemical formation and destruction rates. The 
total OH + formation rate is 



ToH+ - ^N,J - "H^OH+f; 



(A.2) 



where ( is the ionization rate and 6oh+ is the efficiency with 
which ionizations are transfered to OH + production. The de- 
struction r ate Ay.j is assumed to be independent of the consid- 
ered level (Brudere r et al.L 1201 Oal) . and for OH + is given by 



AjV,/ = A H+ = «HlYe£0H+|e + 0.5/h 2 &OH+|H 2 )- 



(A.3) 



In steady state, Aoh+ = roH+/("HAfOH+X which is expected to 
be < 10~ 5 s~\ Since this value is much lower than the AnjiN'J' 
values (> 0.05 s _1 , Table [TJ, the OH + molecules are destroyed 
once they are in the ground (N, J) = (0, 1) level, and eq. (IA.U 
yields for excited N > levels 



N,J;N-\,J' 



Tjv,j + n N+ ij>A 



N+l,J';N,J- 



For the ground (/V, J) = (0, 1) level we can write: 



«0,lA O H 



= r ,i + ^n ur A l 



(A.4) 



(A.5) 



Equations (lA.4t and (lA.5b are easily solved iteratively by start- 
ing from a high TV-level down to the N - state. Results show 



that for^oH+ > 10^, f < 10~ 12 s -1 , and 6 oh+ ~ 0.5, the three 
/V = 1 — > components have an excitation temperature T ex < 7 
K, implying that less than 0.3% of OH + molecules are in excited 
levels. Therefore, excitation by dust emission is much more ef- 
ficient than formation pumping in these sources. A similar anal- 
ysis was done for H20 + , with also similar results: T ex < 7 K 
for the ground state transitions, and less than 1.5% of H20 + 
molecules in excited levels. 

The flux predicted in the OH + /V — > TV - 1 transition as 
derived from formation pumping alone can also be easily esti- 
mated. By summing over all /-levels in eq. (1A.4I >. one obtains 



N'>N J 



(A.6) 



thus the number of N — > A' - 1 line photons emitted per unit 
volume and time is equal in steady state to the accumulated 
formation rate in all N' > N levels, and is independent of 
the A„;-values. Writing the volume of the emitting region as 
V = nR 2 Nulnn, where A^h and «h are the column and density of 
H nuclei, respectively, and R is the effective radius of the region, 
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the flux in the N — > N ■ 
pumping is given by 



1 fine-structure lines due to formation 



hv 



N^N- 1 



N^>N- 1 



R 2 N H v 



4n H D 2 



Yr 



N'J, 



(A.7) 



where D is the distance to the source. The double sum in 
eq. ( IA.7t is evidently lower than the total OH + formation rate 
(roH+X an d tnen 



hv N -, N -iR 2 N H roH+ 



4n H D 2 
and using eq. (1A.2I >. 

hv N ^ N ^R 2 NHeoH+( 



(A.8) 



AD 2 
1.3 x 10" 23 e OH+ 

/ N H 

Ux 10 23 cm- 2 



150yum\/ R 



)(: 



50 pc 

) 



72Mpc 
D 



Wcm . 



(A.9) 



The reference value is ~ 500 times lower than the absorbing flux 
of the N = 2 <— 1 transition in Arp 220 (Table[T]i, indicating that 
formation pumping has a negligible effect on the observed line 
fluxes. For NGC 4418, the upper limit is 3 x 10~ 23 WcirT 2 for 
R = 20 pc, also much lower than observed. 

Appendix B: Excitation of H 3 + through formation 
pumping 

As discussed in Sect. 13. II the rates for H30 + excitation through 
collisions are unknown and therefore the excitation mechanism 
of the high metastable levels cannot currently be determined. 
Nevertheless, here we attempt to constrain under what condi- 
tions formation pumping can account for the observed excita- 
tion. By assuming a range of reasonable collisional excitation 
rates and some simplifying but illustrative assumptions, we show 
here that formation pumping can play an important role in pop- 
ulating the highly excited metastable levels from which absorp- 
tion is observed in Arp 220. 

We assume that the formation distribution can be described 
by a format ion temperature, 7f orm , so that the formation rate T, 
in level i is (Bruderer et al.L 1201 0a): 



r, = r H3 o+ x 



gjexp{-Ej/T foim } 

G(7W) 



(B.l) 



where Fh,o + is the total H30 + formation rate, g, and £, are 
the degeneracy and energy of level i, respectively, and Q is 
the partition function. H30 + is expected to be formed through 
H 2 + + H 2 -> H 3 + + H, with an exothermicity of 1.7 eV. 
Imposing equal momentum of the products in the rest frame, 
and assuming energy equipartition of H30 + between rotational, 
vibrational, and translational energy, we roughly estimate that 
~ 0.08 eV per formed H30 + molecule will go into rotational 
energy, implying Tf orm ~ 800 K. We adopt in the following 
7W = 1000K. 

After a molecule is formed in a J P K level, it will cascade down 
very quickly to the metastable one, where it will remain un- 
til a collisional event induces a transition to another level, or 
until the molecule is destroyed through recombination with an 
electron. If the latter process dominates, the population in the 
metastable levels will be determined by the formation process. 
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Fig. B.l. (a) The H30 + rotation diagram (triangles), calculated 
in the optically thin limit. The colored squares (red: ortho, blue: 
para) show the expected column density distribution for a for- 
mation temperature of Tf orm = 1000 K, and neglecting colli- 
sional relaxation. The two solid lines correspond to rotational 
temperatures of T rot = 250 and 600 K. (b) The rotational tem- 
perature of the 10| metastable level of H30 + relative to the 
5 J one as a function of the electron abundance ^(e) for two 
adopted values of T gas (50 and 200 K) and four assumed val- 
ues of Jfcgjp = 3 x 10" 11 (curves 1), 1 x 10 -10 (2), 3 x 10 -10 (3), 
and 1 x lfr 9 cm 3 s" 1 (4). 



However, the rotational temperature T rot of the metastable levels 
will be lower than 7f olm , because the former is determined by 
the J -accumulated formation rates in the different /f-ladders. 
In Fig. IB. lb . the colored symbols show the expected popula- 
tion diagram obtained by ignoring collisional relaxation. While 
there is not a unique T rot , there is a tendency for it to increase 
with Ei, and the expected distribution can be roughly fitted with 
T m ~ 250 and ~ 600 K for E, < 300 K and E, > 300 K, 
respectively. Th is resembles the double distribution found by 
iLis et alj d2012l) toward Sgr B2. The triangles in Fig. IB. lb show 
the columns obtained toward Arp 220 by assuming that the ob- 
served lines are optically thin (Sect. l23T >. 

In order to check under which conditions the formation pro- 
cess can generate the observed high excitation of the metastable 
H30 + levels in environments with moderate 7^, we have made 
calculations for the equilibrium populations uk of the metastable 
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levels by solving the rate equations 
dviK 

= i k + «H 

K' 



dn K v / col 



- n K n n 



V K' 



+ |e 



= 0, 



(B.2) 



where Yk denotes the accumulated formation rate in all lev- 
els of the ^-ladder. The rate coefficient for H30 + destruction 
through recombinatio n is taken from the UMIST06 database 
(IWoodallet all 120071) : 



%,o+|e = 4.32 x 10~ 7 x 



T 

gas 

300 K 



-1/2 



3 -1 

cm s . 



(B.3) 



We treat the collisional relaxation in a rather simplistic, but still 
illustrative way: only collisions between a given metastable level 
and the neighboring upper and lower (ortho or para) metastable 
levels are considered, with an effective rate of collisional de- 
excitation from level K to level K' given by «h k c £ l K , , assumed 
independent of K. The excitation rate is fixed by the requirement 
of detailed balance at the adopted T gSLS . 

We solved eqs. ( lB.2t for all metastable levels by assuming an 
H30 + ortho-to-para abundance ratio of 1, and for four adopted 
values of k°£ K ,: 3 x 10 -11 , 1 x 10~ 10 , 3 x 10" 10 , and 1 x 10 -9 
cm 3 s _1 . In steady state, Yk «h and under our simple as- 
sumptions the r rot values are independent of «h- As an exam- 
ple of this approach, Fig. IB. lb shows the T mt value between the 
K - 10 and K = 5 metastable levels as a function of Xe- We 
adopt here T gs , s = 50 and 200 K, much lower than 500 K. Since 
both the recombination rate increases, and k°° l K , is expected to 
decrease with decreasing r gas , we may expect the effect of for- 
mation pumping and T Iot to be the highest for the lowest r gas . 
For r gas = 50 K, TVot is much higher than r gas even for relatively 
low x s ~ 10~ 6 . r ro , ~ 500 K is then obtained for ^ e ~ 10~ 4 , 
typical of translucent clouds, provided that k^^, < 10 -10 cm 3 

s _1 . Higher r gas = 200 K implies a lower recombination rate 
and a presumably higher k c £ K ,, with the consequence that J rot 
decreases for the same^f e . From the excitation point of view, we 
thus roughly expect two opposite solutions to explain T rot > 500 
K as observed in Arp 220: low r gas (< 100 K) and, of course, 
high T gas > 500 K. 



o 
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Fig. C.l. Predicted OH + abundance (upper panel), OH + /H 2 + 
and OH + /H30 + ratios (mid panels), and H + and e~ abundances 
(lower panels), as a function of the molecular fraction /h, (up 
to approximately /™ x , see eq. d2j). Each curve corresponds to 

a value of £/«h, indicated in the upper panel in units of cm 3 
s~'. Solid and dashed curves correspond to T gas = 150 and 550 
K, respectively. The inferred values of xoh+, Ri, and R2, in the 
nuclear regions of NGC 4418 and Arp 220, and in the C ex tended 
of Arp 220, are indicated with the shaded dark and light regions, 
respectively. In the lower panel, squares and triangles show the 
predicted xh+ and^ e . 



Appendix C: The formation and abundance of OH + 

Here we present the full set of calculated abundances of OH + , 
H20 + , and H30 + , together with the corresponding abundances 
of electrons and H + , based on the methods discussed in Sect.|4] 
Some analytical expressions for OH + production rates are com- 
pared, showing that the most likely path for OH + formation is 
+ + Hi, with the contribution from the charge exchange reac- 
tion H + + OH and from direct OH ionization being of secondary 
importance even in the nuclear regions, where high OH abun- 
dances are inferred (Paper I). We also show that H + is most 
likely formed through X/cosmic ray ionization, rather than via 
the C + + OH -> CO + + H -> CO + H + FUV/chemical path. 

The predicted OH + abundance relative to H nuclei, xoh* , and 
both Ri and R2 are shown in Fig. IC.ll as a function of /h 2 for 
several values of £/«h, together with the values inferred in the 
nuclear regions of NGC 441 8 and Arp 220 (shaded dark regions) 
and in the C ex tended of Arp 220 (shaded light regions). The abun- 
dances of H + and e~ are also shown in the lower panel. 



The values relevant to the nuclear regions, xoh+ — (1 - 3) x 
10~ 8 , R] ~ 1 - 2.5, and R2 ~ 1, can be obtained for very dif- 
ferent values of £/«h, provided that /h, attains the appropriate 
regime. For example, values compatible with observations are 
found for <T/« H ~5x 10~ 19 cm 3 s _1 , /h 2 ~ 0.1, and high r gas . 
For low values of £/«h, however, the range of fa, compatible 
with observations is very narrow. This range broadens for higher 
£/«h and fn 2 , which we favor owing to the inferred high column 
densities. 

According to the equilibrium xn* and x^ we can compare 
the OH + volume formation rates due to Q + + Ht, Toh +^+ih^ and 
due to O + H3 , roH+;H+|o> at high fn 2 (see iLisztl 12007, for diffuse 
clouds). In what follows, kx\Y is the rate coefficient for the X + Y 
reaction. On the one hand, roH+ ; o + |H 2 = 0.5/h,m^o + ^o + |h 2 > 
where xo* ~ ^H+|o^o»^H+(0.5/H 2 /t o +|H 2 +(1 -/h 2 )^o+|h) _1 . Using 
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Table C.l. Relevant reaction rates involved in the OH + , H20 + , 
and H30 + formation and dest ruction. Rate coefficie nts are taken 
from the UMIST06 database (IWoodall et aUl2007l) . 



Reaction 



Rate coefficient 
(cm 3 s~') 



H + + O -> + + H 


^H+IO = 


7.31 x 10- 10 7/ 3 ° 2 3 exp{-225.9/7:) 


+ + H -> O + H + 


&0+IH = 


5.66 x 10- 10 r 3 ^ 36 expj8.6/7;} 


+ + H 2 -> OH + + H 


feo+|H 2 = 


= 1.7 X 10- 9 


OH + + H 2 -» H 2 + + H 


&OH+|H 2 


= 1.01 x 10~ 9 


H 2 + + H 2 -> H 3 + + H 


£h 2 0+|H: 


, = 6.4 x 10- 10 


H + + OH -> OH + + H 


^H+|OH : 


= 2.1 x 10-" 


H + + H 2 -> H 2 + + H 


&H+|H 2 


= 6.9 x 10~ 9 


H+ + H 2 -> H+ + H 


^Ht|H 2 " 


= 2.08 x 10~ 9 


H+ + -> OH + + H 2 


Woh+ =8.4x 10" 10 


H+ + -> H 2 + + H 


^H+|0:H 2 


0+ = 3.6 x 10- 10 


H+ + OH -> H 2 + + H 2 


^H+|OH ; 


= 1.3 x 10- 9 


H+ + H 2 -> H 3 + + H 2 


^H+|H 2 


= 5.9 x 10-" 


OH + + e -> O + H 


KOH+|e = 


= 3.75 x 10- 8 r 3 -°j 5 


H 2 + + e — > products 


^H 2 0+|e 


= 4.3 x 10- 7 T 3 -f 


H 3 + + e — > products 


«H 3 0+|e 


= 4.3xio- 7 r 3 -» 5 


H 3 + e — > products 


&H+|e = 


6.7 x io- 8 r 3 -°j ?2 



OH+;0+|H 2 



2.5 X 10 _1 

Xo° 



X 



XH+ 



3 x 10" 4 6x10 



1+ 0.55(1 -/ H2 )//h 2 

-3 -1 

-cm s 



(CI) 



On the other hand, r OH +;H+|o = h^To^hjiq^hj, where xn+ 



k H l le (fH 2 /Xe)(£/ n H)- Therefore 
r H + ;Htlo * 1.7 XlO" 19 



X «y X 



Xo° 



X 



3 x 10- 4 



X 



-3 -1 

cm ~s 



3.3 x 10 3 " 2xl0- 17 ^„ • 
Taking the limit fij 2 w 1 in eq. dC.U . one obtains 

ToH+iO+IH, . . Xe/fa 2 XW 2 X 10" 17 



OH+;0+[H 2 
"0H + ;H+|0 



15x 



x ■ 



3 x 10- 4 6x10 



■x- 



cm 3 s 1 



?/n H 



.(C.3) 



The reference values in eqs. ( IC.U - dC.3b are based on £/«h 
2 x lO -17 cm 3 s~' at T g = 150 K, and indicate r H+ ; o+|H 2 » 
Toh+;H+|o; increasing £/hh would increase both^-e+ and;^ e , thus 
further increasing the above ratio. The O + H3 mechanism only 

becomes competitive for low ^/«h 5xl0 _19 cm 3 s _I and/n 2 ~ 
1, but the resulting OH + abundance is then much lower than the 
inferre d one (F ig.l C.U a nd the OH + /H2Q + ratio becomes lower 
than 1 dHollenbach et al.Ll2012l) . 

In the nuclear regions of both NGC 4418 and Arp 220, high 
reservoirs of OH and H2O are found, with estimated abundances 
Xh 2 o ~ (0.5 - 1) x 10- 5 and OH/H 2 ~ 0.5 (Paper I). The OH + 
volume formation rate due to H + + OH is given by roH+;H+|OH = 
"h^h+A'oh^h+ioh, and so 



OH+;H+tOH 



3.8 x 10~ 19 4 x 



Xou 
3 x 10" 



Xh* 



6x10- 



■ cm~ 3 s~{C.4) 



which is significantly lower than roH+;0+|H 2 in eq- dC.ll ) if there 
is enough free oxygen. However, a significant fraction of O in 



a 7/300 is 7//300 K; T denotes r gas . 

the UMIST06 reactions rates given in Table|CJjat j g = 150 K, 

1 



these high density regions is expected to be locked into CO, 
H2O, OH and O2. Depletion of free O into molecules in re- 
gions dominated by hot core chemistry (Paper I) would thus 
reduce the OH + formation rate in eq. dC.ll l. thus increasing 
the potential relative importance of roH+;H+|OH- There is, how- 
ever, an observational drawback for the potential importance 
of H + + OH. The rate of H z O + formation due to H + + H z O is 
higher than roH+;H+|OH, because both^HjO > Xou (Paper I) and 
£h+]h 2 o > &h+|oh (Table IC.U . Thus, if the OH + were formed 
via H + + OH, one would also expect an additional, important 
contribution to H20 + formation due to H + + H2O, and then Ri 
would fall below the inferred value of 1 - 2.5. Likewise, direct 
ionization of OH due to UV photons is expected to be of sec- 
ondary importance, because H2O has a higher cross-section for 
ionization. As we argue in 34.51 the observed lines of molecular 
ions are most probably formed in regions where the formation 
rate is enhanced due to the presence of high quantities of atomic 
oxygen, and where the ions survive sufficiently long due to the 
presence of moderate densities. 

According to the above discussion, model results for xoh* 
and Ri relevant for the current observations can be basically un- 
derstood in terms of the H + -> + -> OH + -> H 2 + -> H 3 + 
sequence, but how is the H + produced? We have been assum- 
ing in the analysis that the H + ionization is caused by X/cosmic 
rays, but the possibility of e nhanced H + production by c hem- 
ical routes, as di s cussed by ISternberg & Dalgarnol d 1995b and 
Hollenbach et a 

D d20ll . is worth considering. According to 
their models for high density PDRs illuminated by strong FUV 
fields, in the hot transition region between the H 1 and H2 zones, a 
major source of H + formation is given by OH + C + — > CO + + H 
followed by CO + + H -> CO + H + . OH would be produced via 
the endothermic O + H2 — » OH + H, that is enhanced provided 
that a fraction of H2 is vibrationally excited. The volume rate of 
H + production is given by r H+; 0H|c+ = n^k 0Hlc + XohXc* ( 1 +a)~\ 
where a stands for the chance that CO + reacts with H2 rather 
than wit h H 1. In the limit a = 0, and using the value for £oh|c+ 
given bv lHollenbach et all (1201 2l) . 



(C.2) rH+ ; oH|c+ - 2 x 10 



8 »uX 



Xc+ 



1.4 x 10" 



Xon 
5 x 10" 



'(C.5) 



for r gas = 300 K; the reference value for ^oh is taken from 
Paper I. This is equivalent to a cosmic ray ionization per unit 
density of £/«h ~ 2 x lO -18 cm 3 s _1 , and predicts a peak 
of OH + abundance of ~ 1Q- 8 as well as R\ and R2 value s 
of ~ 1 (ISternberg & Dalgarnol Il995t iHollenbach etalL 120121) . 
as roughly inferred in both NGC 4418 and Arp 220. This 
FUV/chemical route for the production of H + has in the present 
case, however, one main drawback: the nuclear OH + is observed 
over a range in visual extinction of AAy ~ 200 mag, while the- 
oretical predictions indicate a peak thickness in the OH + abun- 
dance of AAy < 1 mag where C + is still abundant and r gas re- 
mains high. It seems very unlikely that FUV radiation maintains 
basically all carbon ionized, and a relatively low H2 fraction, 
over columns of AAy ~ 200 mag along the line of sight and 
over spatial scales of 20 - 30 and 80 - 150 pc in NGC 441 8 
and Arp 220, respectively (see also models bv lAbel et al.l l2009). 
Even if 10 6 O stars were packed into a region of ~ 10 pc ra- 
dius (NGC 4418), the column density between neighboring stars 
would be ~ 10 23 cuT 2 , and the quickly absorbed FUV radiation 
would only affect relatively compact regions around the stars. 
In summary, the inferred columns of O-bearing cations favor a 
source of ionization that can penetrate through large columns of 
gas, indicating the importance of X and cosmic rays. 
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The dependence of xoh+ on both /h, and £/«h in Fig. IC.ll 
can be understood as follows. For given /h 2 , ^oh+ scales lin- 
early with £/«h for relatively low values of £/«h, but saturates 
for sufficiently high values of f /«h- Saturation occurs when the 
destruction of OFF is dominated by recombination, rather than 
by reaction with H2 (producing H20 + ). In the saturated regime, 
further increase of £/«h increases by the same factor the OH + 
formation rate and the OFF destruction rate per unit of OFF 
molecule, leaving xoh+ unchanged. The condition for saturation 
is given by ^-e^OH+|e » 0.5/h 2 &oh+|h 2 , where k OH * ]e ~5x 1(T 8 
cm 3 s -1 at 150 K, and &oh+|Hi = 10~ 9 cm 3 s _1 . The above con- 
dition yields Xe 2> 10~ 2 x fu 2 . Once saturated, xon+ can t> e a P _ 
proximated analytically by 

sat „ k W\Q 0.5/h 2 ^ O+ |H 2 

° H+ ~ £ H+|e 0.5/H 2 ^O +t H 2 +(l-/H 2 )^O + |H 

* 7.9 x 10~ 7 x x — , (C.6) 

1+ 0.55(1 -/ H2 )//h 2 3x10- 4 

giving, for T gas = 150 K,^ + « 2.7 x 10~ 8 and 3.4 x 10~ 7 for 
/hj = 0.02 and 0.3, respectively, in rough agreement with the 
values in Fig. IC.ll 

The dependence of xon+ on fa 2 is different in the saturated 
and unsaturated regimes. According to eq. iC.6i . yrm+ increases 
with /h 2 in the saturated regime, because the + + H2 reac- 
tion producing OFF is increasingly favored relative to + + H, 
which transfers the charge back to atomic hydrogen. However, 
in the unsaturated regime xoh* decreases with increasing /h 2 
(Fig. IC.U . In this case, OH + destruction is dominated by 
OFF + H2 and^0H+ is thus analitically aproximated as 

^unsat ^ fc H+|Q ^ fcp+|H 2 *H+ ^ ^ 

^OH+|H 2 0.5/h 2 ^O+|H 2 + (1 - /h 2 )^o+|h 

« 1.5 x 10- 8 x l - x x ^ (C.7) 

1+0.82/h 2 3x10- 4 10- 4 

Therefore, x^aa* decreases with increasing /h, . For given £/«h 
(i.e. for each curve in Fig. IC.U , a maximum of OFF is reached 
between the saturated regime (low /h 2 ) and the unsaturated one 
(Mgh/ H2 ). 
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